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The backscattering of fast (0.3, I, .'$, 6, 9, and 14 Mev) neutrons from a semi-infinite 
water medium has been calculated by the Monte Carlo method. The information obtained 

includes the joint angular and spectral distribution of the reflected neutrons, the dependence 
of the number albedo and energy albedo on the source energy and obliquity, and the con- 
tributions to the albedo of successive orders of scattering. The spectra were calculated 
down to epithermal energies (~0.5 ev). The results for each case are based on the analysis 
of 3,000 neutron histories, generated with the use of an EBM-704 computer. In the random 
sampling, elastic scattering from hydrogen and oxygen, inelastic scattering from oxygen, 
and absorption due to n-ct and n-p processes were taken into account. The cross sections 
for some of these interactions are not well known. Parallel calculations with different 
assumptions about the cross sections were? made in order to estimate how sensitively the 
albedo depends on the cross sections. The paper includes a self-contained description of 
the Monte Carlo method, its application to the calculation of radiation diffusion and in 
particular to the neutron albedo problem. Emphasis is placed on the technique of cor- 
related sampling which makes possible an accurate estimate of albedo differences resulting 
from different assumptions about the cross sections. The random sampling computations 
were supplemented by analytical calculations of the single-scattering albedo. This was 
useful for an understanding of the Monte Carlo results because a considerable fraction of 
the reflected neutrons return after only one collision. 

1. Introduction 
1.1. Statement of the Problem 

Neutrons that are incident on an extended medium make one or more collisions with 
nuclei in each of which they are deflected and lose a fraction of their energy. Tims they 
perform a "random walk" until they arc either absorbed or re-emerge front the medium with 
reduced energy. The probability that an individual neutron will eventually re-emerge from a 
semi-infinite medium (i.e., be reflected) is called the albedo, or the differential albedo if the 
energy and or direction of the reflected neutron are specified. 

The purpose of this paper is threefold: (1) To calculate the differential albedo and related 
quantities for neutrons incident on a semi-infinite water medium; (2) to give a self-contained 
exposition of the methods of calculation used (Monte Carlo, single-scattering analysis); (3) 
to analyze how the albedo depends on various factors, particularly the cross sections for elastic 
and inelastic scattering and absorption. 

1.2. Background 

Exact calculations of the neutron albedo have so far been limited to situations in which 
the slowing down of neutrons could be disregarded. With this limitation, and with further 
restrictive assumptions, Halpern, Lueneburg, and Clark [1] 2 have given an analytic solution 
bused on transport theory. A more general solution of the one-velocity problem is also avail- 
able from Chandrasekhar's [2] theory of radiative transfer. Fermi [3] and more recently 
Grosjean [4] have given approximate solutions which are in good agreement with the more 
rigorous calculations. 

When the complicated energy-dependence and directional dependence of the collision 
processes must be taken into account, an exact analytical treatment becomes all but impossible. 



1 Supported jointly by Armed Forces Special Weapons Project, OCDM, and the Department of the Navy (Bureau of Yards and Docks). 

2 Figures in brackets indicate the literature references at the end of this paper. 
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Spinney [5] has given an approximate treatment based on age theory which is applicable 
mainly to energies lower than those considered in the present work. 3 Another possible approach 
would be to incorporate the exact one-velocity treatment of reference [2] into an energy-group 
scheme similar to the multigroup diffusion schemes used in reactor calculations. This has not 
been attempted so far. 

The Monte Carlo method recommends itself as the simplest and most promising approach. 
It has the advantage that the complexity of the calculations is not essentially increased by 
adding to the number of variables or by introducing complicated cross sections. Increased 
detail merely creates a data-handling problem which requires a computer with sufficient ca- 
pacity. A disadvantage of the Monte Carlo method is its purely numerical character which 
makes necessary a very careful scrutiny of the results before one can understand their signifi- 
cance. We have found that such an analysis is greatly facilitated by parallel analytical 
calculations of single scattering. In spite of its promise, the Monte Carlo method has not 
yet been used extensively for albedo calculations. Some results for the backscattering of 
0.89-Mev neutrons from water have been given by Foderaro and Obenshain [6], but their 
work was mainly concerned with transmission through slabs. 

One reason for the scarcity of Monte Carlo applications to the albedo problem appears to 
be the lack of adequate information on neutron collision cross sections, which may have dis- 
couraged extensive numerical work. 4 For few if any elements (other than hydrogen) are the 
cross sections known in the detail and with the accuracy desirable. Only the knowledge of 
the total collision cross section is in definitive form for most elements. Information about 
the contributions of different processes to the total cross section is reasonably adequate. Much 
less is known about differential cross sections. The angular distribution of elastic scattering 
deflections — which is highly anisotropic at high energies — is not well known for many elements. 
Even less is known about inelastic scattering. There is uncertainty as to the nuclear levels 
which are excited, the relative excitation probabilities, and the angular distribution of the 
inelastic scattering deflections. 

We have chosen water as the medium in our calculations because the cross sections for 
hydrogen and oxygen are known better than those of most other elements. In addition to 
being of interest in their own right, the results for water will also be typical of those expected 
with other hydrogeneous substances. Thus the present work should provide guidance for 
future calculations for media such as concrete which are held up now by lack of cross section 
data. 

2. Variables of the Problem 

We consider the reflection of neutrons from a semi-infinite medium occupying the region 
2>0. Figure 1 shows a typical trajectory of a neutron which is incident on the medium with 
obliquity 6 and which is reflected with obliquity 6. These obliquities are angles between the 
direction of motion of the neutron and the positive 2-axis. Thus 6 lies between and tt/2, 
and 6 lies between ir/2 and ir. 

We denote by A(E,d;E ,d ) sinddddE the probability that a neutron incident with energy 
E and obliquity o will be reflected from the semi-infinite medium, after one or more scatterings, 
with an energy between E and E-\-dE and an obliquity between 6 and d+dd. The function 
A(E,d;E ,S ) is called the differential albedo. Alternatively, instead of regarding the albedo as 
a probability, we may consider it to be the ratio of the average number of neutrons reflected 
to the number of neutrons incident. Thus if a beam of N neutrons is incident with energy 
E and obliquity d , and if N(E,0) sinddddE neutrons are reflected on the average, then 



3 See, however, sec. 7.8. 

4 For a survey of the neutron cross section situation in relation to transport calculations, see the recent book by Goldstein [7]. 
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By integrating the differential albedo one obtains the number albedo A N (ratio of the number 
of all reflected neutrons to the number of incident neutrons): 

A N (EM= [ H (IkV smedOA(E,e;E ,e ). (2.2) 

JO J 7T/2 

Similarly, the energy albedo A E (ratio of the average reflected to the incident energy) is defined as 

A S (E ,6 ) = I °dE ( * miM%A(E,0-,E o ,eo). (2.3) 

JO Jtt/2 J^O 

We shall have occasion to consider the probabilities of reflection after exactly 1,2,. . . n 
collisions. The corresponding albedos will be indicated as A 0l) (E, 6;E o ,0 o ), A N in) (E O) O ) and 
A E in) (Eo,0 ),n=lj2, .... The case n=l (single-scattering albedo) is particulary important. 




Figure 1. Neutron trajectory. 



3. Neutron Cross Sections of Hydrogen and Oxygen 

3.1. Integral Cross Sections 

The most important parameter characterizing the neutron transport process is the mean 
free path between collisions, \(E). It is inversely proportional to the number ~^of target 
nuclei per unit volume and the total collision cross section <r(E), 



\{E) = 



1 
Jfo(E) ' 



(3.1) 



The reciprocal of the mean-free path (the attenuation coefficient) is additive for a homogeneous 
mixture of different nuclei. For water, 



^ = ~+^=Jk<ru(E)+^ 0x o- 0x (E) y 

AHoO A H Aqx 



(3.2) 



where ^=6.69X10 22 cm- 3 and^ x =3.34X10 22 cm- 3 . 

We have taken the values of <r H (E) and a 0x (E) from the standard Brookhaven compilation 
[8]. The attenuation coefficients for hydrogen and oxygen thus obtained are shown separately in 
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figure 2 as functions of the energy. The plots are made in histogram form corresponding to the 
approximate manner in which the data were used in the machine computations. It can be seen 
that the energy-dependence for hydrogen is smooth, but that the histogram for ox}^gen has 
several sharp peaks and irregularities at high energies. 

Next we consider the relative contributions of various processes to the total collision cross 
section in the energy region of interest (0.5 ev to 14.0 Mev). For hydrogen, only elastic scat- 
tering is important in this region. The (n,y) cross section is smaller than 0.1 percent of the total 
cross section at 0.5 ev, and is even less significant at higher energies. 
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Figure 2. Attenuation coefficients for hydrogen and oxygen. 
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For oxygen, the processes that mallei- are clastic scattering, inelastic scattering (above 
— 6.4 Mev), and absorption due to (n,a) reactions (above«4 Mev) and (n,p) reactions (above — 
10 Mev). Cross sections for these processes, derived from the best experimental evidence and 
theoretical considerations, have been compiled by Lustig, Calos, and Goldstein [9]. 5 In table 
3.1 we list the values of the various cross sections at selected energies to indicate their relative 
importance. For use in machine calculations, we have fitted cross-section ratios with poly- 
nomials in the energy variable. The representations are a smoothed-oul version of the data of 
reference [9] accurate to 5 percent or better, which is sufficient considering the uncertainty in 
the data. 

Table 3.1. Neutron cross serf ion* (barns) for oxygen (from referenci [9]) 



Energy 
Mev 


ff o>. 

Total 


Vox 

II, a 


%,0x 

II - P 


<V,o* 
Inelastic 


11 
(i 


1 . 58 
1.08 
1.48 
1.30 


0. 31 

'. 19 


0.09 


0. 5 ) 
.21 


;< 















Ratio of inelastic to total collision cross section: 

-1 .05173+0.24196£ r — 0.013666£ ,2 +0.00024009£ ,3 (6.4 Mev <./<;< is Mev) 



0Ox 



Ratio of absorption to total collision cross section: 
^o'+^.o' ^— 0.40433+0.13254g— 0.0085614g 2 +0.00017413g 8 (4 Mev< E< 18 Mev). 



Y (3.3) 



°"<)x 



The probabilities that a given collision will be with a hydrogen or an oxygen nucleus are 

X H (E)/\ H (E) and X H (E)/\ x(E), respectively. A collision can have four possible outcomes: 
1 . Wit h probability 

Xh 2 o 



J>r- 



(3.4) 



the neutron is scattered elastically by a hydrogen nucleus. 
2. With probability 



Vi- 



Xh 2 o 



°"o.\ °V, Ox °"a,Ox °"/>.(>\ 

— ; 
<?Ox 



the neutron is scattered elastically by an oxygen nucleus. 



3. With probability 



lh= 



^H 9 a,,, 



2° OV.Ox 



X() X 0- ()x 



the neutron is scattered inelastically by an oxygen nucleus. 
4. With probability 

^ii.,<> o- ( , )0x +<rp )0x 

A ox »Ox 

the neutron is absorbed. Note that ^1+^2+^3+^4= 1. 



(3.5) 



(3.6) 



(3.7) 



5 Duo to theoretical and experimental uncertainties, the cross sections given in this publicai ion are subject to revision. We understand that 
Dr. Goldstein and his associates arc now engaged in making an improved evaluation. 
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3.2. Differential Cross Sections 

The outcome of a collision in which a neutron of energy E is scattered can be specified in 
terms of the probability distribution jj(Q)E) of the scattering deflections 9 C in the center-of- 
mass coordinate system. In conformity with the notation in the preceding subsection, we 
shall use indices j=l, 2, and 3 to indicate the distributions for elastic scattering from hydrogen, 
elastic scattering from oxygen, and inelastic scattering from oxygen, respectively. 

For elastic scattering from hydrogen, the distribution is assumed to be isotropic, 6 i.e., 

/i(G c ;S) = l. (3.8) 

The angular distribution for elastic scattering from oxygen is also isotropic at energies 
below about 300 kev, but markedly anisotropic at higher energies, particularly in resonance 
regions. Above 300 kev we have used the compilation of reference [9] which gives the coef- 
ficients / 2 ,zCE') of the Legendre expansion 

/ 2 (e c ;£)=S ^fi-fUE)P l (cos C ). (3.9) 

l * 

(Note that the normalization is such that / 2 , (£') = 1.) Typical angular distributions at various 
energies are shown in figure 3. 

The angular distribution for inelastic scattering from oxygen depends not only on the 
energy of the neutron before the collision but also on the excitation energy which the struck 
nucleus receives. In O 16 , the four lowest levels are 6.06, 6.14, 6.91, and 7.12 Mev above the 
ground state, the next is at 8.6 Mev and there are numerous levels between 9 and 14 Mev [11]. 
It appears from the discussion in reference [9] that neither the relative probabilities for the 
excitation of these levels nor the corresponding angular distributions are known. In order 
to explore the effect of inelastic scattering on the albedo, we have done three sets of calculations 
on the basis of the following assumptions: 

(1) The nucleus is excited to a level 6.1 Mev above the ground state. The angular dis- 
tribution is isotropic; i.e., 

/ 8 (e c ;E) = l. (3.10) 

(2) The neutron loses all its energy in an inelastic scattering, so that it is, in effect, 
absorbed. 

(3) An inelastic scattering is treated as if it were an elastic scattering; i.e., the excitation 
energy of the nucleus is set equal to zero. 

The assumptions of case (1) seem to us to be the most realistic, and may be expected to 
yield reasonable answers when the neutron source energy is not too much above 9 Mev. The 
assumptions of case (2) lead to an underestimate of the albedo, and the assumptions of case 
(3) lead to an overestimate of the energy of the reflected neutrons. 

4. Single-Scattering Albedo 

4.1. Formulation 

The probability that a neutron, incident on a semi-infinite medium with energy E and 
obliquity 6 , will be reflected with energy E and obliquity 6 after exactly one collision is 

f 00 d? 3 

A^(E,6;E ,d )=\ e-'^-^PjHjiEfiiEoMe- 8 '* iiir/2<6<T } 

Jo *0j=l ( 4 ^ 

= ifO<0<7r/2. 



e This assumption is open to doubt at energies above, say, 5 Mev. See Goldstein, Zweifel and Foster [10]. 
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Figure 3. Angular distribution of elastic scattering from oxygen (according to ref. [9]). 

In this expression, the factor e" s/x o — represents the probability that the incident neutron will 

travel a pathlength s and then undergo its first collision between s and s+ds, \ =\(E ) being 
the mean-free path. The collision can be of type 1 (elastic scattering from hydrogen), type 2 
(elastic scattering from oxygen), or type 3 (inelastic scattering from oxygen). The fourth 
type (absorption) need not be considered. The factor jp } (defined by (3.4, 3.5, and 3.6)) 
represents the probability that the collision will be of type j. The function Hj{E,B)EM is 
called the scattering function and denotes the probability that in a collision of type j the energy 
of the neutron will change from E to E and the obliquity from d to 6. Finally, exp(— s'/X) 
is the probability that the neutron will travel from the point of collision to the boundary 
without making any more collisions (X is the mean free path at the energy E). The pathlength 
to the boundary is 

s' = -s C -^. (4.2) 

cos v ' 
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Using this relation, and carrying out the integration with respect to s in (4.1), we obtain the 
equation 

;=1 COS V — {\qIA) COS (7o 

(4.3) 
=0 ifO<0<7r/2. 

4.2. Scattering Function 

The function Hj(E,d;E ,d ) to be inserted into (4.3) can be derived from the corresponding 
distribution fj(Q c ;E ) for the deflections in the center-of-mass coordinate system. To establish 
the connection, we need the relations between 9 C , the corresponding deflection angle in the 
laboratory coordinate system, and the ratio of energies before and after the collision, E /E. 
These relations are a consequence of the conservation of momentum and energy in the collision 
of a neutron with a nucleus. We state them for the case of inelastic scattering, characterized 
by the transfer of an amount Q i of energy to the target nucleus as excitation. 7 The same 
relations also apply to elastic scattering if we set Qj=0. 

The energy of the neutron after the collision is given by 

and the deflection in the laboratory system is 

/ff, l+M,g,co9e. _ M,+l [E MfKf-1 1% ( 

cose-y-^ 1+M , ~2~ ^E 2 (M,+ 1) V^' ( } 

where M^ is the mass of the target nucleus and where 



*rfW- 



(4.6) 



Let Fj(E;E )dE denote the probability that the energy of the neutron after the collision 
lies between E and E+dE. Using (4.4) we have 

Fj(E;E ) =/ i (e c ;S ) sin e c ^ ^M&Eo fj{Q '' En) ■ (4 " 7) 

According to (4.5), the deflection G is determined by the value of E/E . The joint distribution 
of E and 6 can, therefore, be expressed formally with the use of a Dirac delta function, as 

F/EirTU / „ Mj+1 \~E , M)K)-\ IEq\ , a v 

F 3 (E;E )8 (cos e— f- V^+2^T) V #/ (4 " 8) 

Let the direction of the neutron before the collision be represented by a unit vector 
(sin0 o cos <p , sin0 o sin<p , cos0 o ), and its direction after the collision by a unit vector (sin0 
cos<p, sin0 sin<p, cos0). Then the cosine of the deflection angle is the inner product of the 
two vectors, 8 

cos6 = cos0 o cos0+sin0 o sin0 cos(^ — <?)• (4.9) 

(The angles <p and <p are azimuths with respect to the x-z plane.) If we insert (4.9) into 
(4.8), introduce the abbreviation 

Mj+l [E MfKj-1 lEo ,. 

^ Sa ^-2-V^~2(M^V^ (4 ' 10) 



7 For a derivation, see, e.g., Bethe [12]. 

8 See also figure 5. 
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and carry out an azimuthal average with respect to <p' = (p {) —(p, we obtain the scattering 
function 



1 f "" 

IIj(K,d;K {) ,d () ) =Fj(E;E ) ~ '/</<5(cos 6 cos 0+sin 0„ sin 6 cos <?'— cos a 3 ) 

_ Fj(E;E ) 



(4.11) 



7r{ cos(o,— -0 O ) — cos# } 1 { cos0— cos^-f O ) }^ 

The square roots in (4.11) must not become imaginary, which requires thai 

cos(aj+d ) <cos0<cos(a:,— O ). (4.12) 

In the case of perpendicular incidence (6 = 0), the region (4.12) shrinks to a single value 
cos a h and the scattering function becomes 

H j (E,d;E ,0)=Fj(E;E )8(cos 0-cos a,). (4.13) 

In the limiting case of grazing incidence (0 O = 9O°), the scattering function becomes 

H/E, 6; Eo, tt/2) = , F 'f o ;E °\ =- (4.14) 

x-ysnr 6— COS" aj 

4.3. Spectra and Angular Distributions 

When using the single-scattering results as a guide for interpreting the Monte Carlo 
results, it is useful to have the differential albedo with one of the variables integrated out; 
i.e., the t nergy spectrum 

A a) (E;E ,do)= ( &medeA (1) (E,e;E ,eo), (4.15) 

and the angular distribution 

A^(e;E 0) do)=jdEA a) (E } d;Eo y do)' (4.16) 

In the evaluation of these expressions the range of integration must include all values 
of 6 for given E (or of E for given 6) for which the condition (4.12) is satisfied. Moreover, 
cos0 must be negative. The resultant formulas are somewhat complex. We give them in 
appendix A, confining ourselves here to some qualitative remarks. 

The chief characteristic of the energy spectrum is that it is confined to limits which are 
wide for the case of hydrogen but very narrow for somewhat heavier nuclei such as oxygen. 
The upper limit of the spectrum follows from the condition that the neutron most be deflected 
sufficiently to acquire a velocity component in the negative ^-direction. This implies, as a 
consequence of (4.12), that E cannot exceed 



w _„ f \~M W* -cos 2 flo+sinfl \* ( 

/W-Ao| Mi +1 "/' (4,1 ° 

The lower limit of the spectrum is determined by the conservation of momentum and energy 
which implies, according to (4.4), that E cannot be smaller than 

«-*<3£r)' ' <" 8 > 
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The allowed energy regions between E max and E min are indicated in figure 4 as functions of the 
obliquity of incidence O . 

The angular distribution cannot be described by any one simple law. For perpendicular 
incidence, the distribution follows approximately a cosine law; i.e., it peaks in directions corre- 
sponding to perpendicular emergence from the medium. As the direction of incidence becomes 
more oblique, this peak tends to disappear. In the limiting case of grazing incidence the angular 
distribution becomes flat or may even peak in directions corresponding to emergence parallel 
to the boundary of the semi-infinite medium. The angular distribution under these circum- 
stances also becomes a sensitive function of the distribution J(Q C ;E ) of scattering deflections. 
For details, see appendix A. 




Figure 4. Possible energies which a neutron may have if it is reflected after a single scattering. 

The allowed energy regions for elastic scattering hold for all Eo, whereas the indicated region for inelastic scattering applies only to 2£o=14 
Mev and 6.1-Mev excitation. 

5. Monte Carlo Method 

5.1. Generalities 

Although the general principles of the Monte Carlo method are well known, 9 we give a 
brief review to make this paper self-contained. The usual assumption is made that a neutron 
undergoing multiple scattering travels in a straight line until it makes a collision with a nucleus, 
is then suddenly deflected, travels again in a straight line till the next collision, and so on. 
The trajectory of a neutron is, therefore, completely described by a specification of the initial 
conditions and of the energies, obliquities, and positions of the neutron immediately after 
successive collisions. In the physical process, these characteristics are chance variables whose 
probability distributions are determined by the neutron cross sections. The Monte Carlo 
method consists of imitating nature by playing a game of chance. With the use of random 
numbers, the parameters of model neutron trajectories are calculated according to the appro- 
priate probability distributions. Such a procedure is called random sampling, and a neutron 
trajectory thus generated is usually called a neutron history. 

9 See, e.g., reference [13] for background and bibliography. 
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By examining a sufficiently large number of neutron histories, one can obtain statistical 
estimates of various transport phenomena. For example, if A 7 ,, histories are generated in each 
of which the neutron starts at the boundary of a semi-infinite medium, and if in TV of these 
histories t he neutron event ually emerges from the medium, t hen the ratio N/N is an estimate of 
the number albedo A N . 

5.2. Random Sampling of Neutron Histories 

If we let E n , n , and z n represent the energy, obliquity, and distance from the boundary 
immediately after the n'th collision, the neutron trajectory or history can be described by the 
array 

E , E u . . . , E ni . . . 

Oq, d u . . . , 6 n) , . . (5. 1) 

z , #i, . . . , z nj . . . 

in which E 0i 6 0j and z =0 indicate the initial conditions. The history ends when the neutron 
is absorbed or reflected from the semi-infinite medium. We have adopted an additional termi- 
nation rule in the present work, ending a history when the energy of the neutron drops below 
0.5 ev. In order to generate a history, one must proceed from collision to collision, computing 
(En+ifln+iyZn+i) from (E n ,6 n ,z n ) recursively by random sampling. 

By "random sampling from a distribution f(x)" we mean, in principle, the following 
procedure. A set of numbers, x u x 2 , . . ., is prepared whose frequency distribution is made 
to approximate /(a?) as closely as possible. One of these numbers is then selected at random 
and is called a random variate from the distribution /(.r). It is, in fact, sufficient to prepare 
only a single set of random variates, from a uniform distribution /(aO = l(0<a;<l). These 
numbers, distributed uniformly between zero and one, are called random numbers. When 
suitably manipulated they can be made to yield random variates from any distribution, as 
shown below. 

In a machine calculation requiring many random numbers it would be awkward to prepare 
and store a large set of random numbers. Instead, it is customary to use so-called pseudorandom 
numbers which are generated systematically but have, nevertheless, sufficient randomness 
for Monte Carlo calculations. We have employed the method of congruent ial multiplication 10 
in order to produce a sequence of pseudorandom numbers 

&, &, • • -, & ■ ■ •• (5-2) 

These numbers are fractions (between zero and one) which are obtained by normalizing a set 
of integers x nj 

Z«=2-**x n . (5.3) 

The integers, in turn, are obtained from the following recursive scheme: 



x =c 
sc n+1 =5*a; n mod2 35 , 



} (5.4) 



where c and k are odd integers. The last equation states that .r„, , is congruent to 5 k x n modulo 
2 35 ; m other words, that x n+l is the remainder obtained upon division of b k x n by 2 35 . In the 
course of our Monte Carlo calculations we have used three pseudorandom number sequences 
of the type (5.2), with k=ll, 13, and 15. This will be discussed more fully in section 5.4. 



111 Sec, (\^., the article by J. Todd in reference [13c]. 
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It can be proved (see app. B) that the sequence (5.2) contains 2 33 different numbers which 
are distributed with uniform density between zero and one. The successive random numbers 
\\>n\ iii the Monte Carlo calculation have been obtained by picking successive numbers in the 
pseudorandom number sequence (5.2). This is equivalent to a random choice, because the 
pseudorandom number sequence closely resembles a true random number sequence in the sense 
that there is no noticeable correlation between successive numbers. This lack of correlation 
has not been proved, but is an empirical fact that can be demonstrated by statistical tests. 
In appendix B, section 9, we present the results of such tests for 237,000 pseudorandom numbers 
actually used in our calculations. 

One procedure for selecting a random variate x from a distribution f(x) by the manipulation 
of random numbers depends on the possibility of finding a random variate y(x) which is dis- 
tributed uniformly between and 1. We le1 

r* co 

y(x)=jj(x')dx', (5.5) 

i.e., y{x) is the cumulative probability distribution associated with /(a;). Note that 

probability {?/<]/' } = probability {x>x(y') } = f Kx')dx'=y', (5.6) 

Jx(yO 

so that y is equally likely to have any value between and 1. We therefore identify y with a 
random number £, and determine a random variate x by inverting the equation 



6= 



--£f(x')dx. (5.7) 



A second method, called the "rejection technique" applies to distributions with a finite 
range, sa % y 0<x<a. Let L be the maximum oif(x). The sampling procedure consists of the 
following steps. First, two random numbers, say £i and £ 2 , are selected. If L% 2 <f(a£i) we 
chose the number afi as a random variate from/(x). If L^^>j{a^), we reject the random num- 
bers £i and £ 2 and repeat the testing process with a new pair of random numbers. The 
justification for this selection procedure is supplied by the observation that 

probability {acceptance of x=a%i) ^probability {l4 2 >f((i%i)} = d£ 2 ——f-' (5.8) 

The average number of random number pairs that must be tested before finding an acceptable 
random variate is given by 

1 



xjj {a ^ d ^ 



=aL. (5.9) 



We shall now describe the actual sampling procedure for generating a neutron history; 
i.e., the rule for going from (E n ,6 n ,z n ) to (E n+1 ,6 n+1 ,z n+1 ). 

a. Location of the Next Collision 

The probability distribution of z n +i is 

e -s n+1 ME n ) = £ e cos*„ A(£ n ), (5.10) 



HE n ) \(E n ) 

where 



(5.11) 



' n+1 ~ cos 9 
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is the pathlength traveled by the neutron between the nth and n-\ 1 s( collision- The distribu- 
tion (5.10) is derived from 1 be assumption thai 1 be spatial dist ribution of target nuclei is random 
and that the probability of a collision perunil pathlength is l/\(E n ). Applying the cumulative- 
probability method, we seled a random number $«, and determine s w _| ( from the equation 

e-WJ r^U=<rWX<*„>. (5.12) 

The pathlength determines the position of the next collision, which is given by 

Z n+1 = 2 n + CO$O n S n+ i = Z n — COSd n \(E n ) log £«. (5. 1 3) 

b. Type cf Collision 

The probabilities Pj(j=l, 2, 3, 4) for the occurrence of each of the four possible types of 
collisions are given by (3.4), (3.5), (3.6), and (3.7). To make a random selection, we divide 
the region between and 1 into four intervals proportional to the probability for each process. 
A random number %p is (hen picked and the type of collision is determined by (he interval into 
which & falls: 

0<^</) { : elastic scattering from hydrogen 

/>i<C&</v! Pi '■ elastic scattering from oxygen 

i>>i+/>2<^/3</vl 7>:H />3 : inelastic scattering from oxygen 

Pi + p2+Pz<b<Pi I p2+Pz+P*=l • absorption. 

c. Energy Loss 

If the ■//-)- 1st collision results in a scattering, the energy of the neutron is changed from E n 
to E n+] according to the single-scattering law (4.7); i.e., the probability distribution of E n +i is 

F t (E n+l ; EJ ^M^k / ' (9c -" +i; En) ' (5,14) 

where 

The choice of the energy loss can he made indirectly through the choice of the center-of-mass 
deflection B c ,n+i- When the distribution of deflections is isotropic (lor collisions of type 1 and 
3), we set 

cos 9 c , n+1 = 2£ 7 -l, (5.16) 

where £ 7 is a random number. For elastic scattering from oxygen, the angular distribution is 
anisotropic and is given by the Legendre expansion (3.9). In this case it is convenient to use 
the rejection technique. Provisionally, a deflection angle is chosen according to (5.10). An 
additional random number, %™ , is then selected and a test is made to determine if 

/ 2 {cos- 1 (2^-l);£'„} >L{E n )$\ (5.17) 

where L(E n ) is the largest possible value oif 2 (Q c ;E n ). If the inequality (5.17) is satisfied, the 
provisional random variate given by (5.16) is accepted; if not, the selection procedure is re- 
peated with two new random numbers. 11 



u Using formula (5.9), noting that the range of the angular distribution (in cos0 c ) is a = 2, and using the differential cross sections as shown in 
figure 3, we find that the average number of random number pairs that must be tested to make a selection has the following values: 6.8 at 14 
Mev, 4.3 at 9 Mev, 2.2 at 6 Mev, 2.6 at 3 Mev, 2.8 at 1 Mev, and 1.03 at 0.3 Mev. The sampling efficiency at higher energies could be improved by 
more elaborate schemes, but we have not considered this worthwhile. 
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d Chang: cf Obliquity 

The prescription for selecting the obliquity after the 7i+lst collision can be obtained in 
two different ways: by a direct geometrical argument, or by an application of the scattering 
function derived in section 4. Starting with the latter, we note that, according to (4.11) the 
probability distribution *of d n+ i (for fixed B c ,n+i and n ) is given by 



where 



7r{cos(o^— 0J — cos0 n+ i}2{cos0 n+1 — cos(a, n +0»)}* 
2 > E n 2(M,+ 1) ME n+1 ' 



cos a Jn = 



and 



2(M, 

cos (a jn + 6 n ) < COS0 W+ ! < cos (a Jn —d n ) . 
Applying the cumulative-probability method, we pick a random number, &, and set 

sinddd 



(5.18) 

(5.19) 
(5.20) 



1-{i= 



1 r<* 1n +e n 



=1 — cos 

TV 



{cOs(a /w — d n ) — COS0}^{cos0— COs(a w + yn )}i 
_! f COS0 n+ i — COSO£y n COS0, 

\ sina^ sin0„ 



■}■ 



The inversion of this equation yields the rule for calculating W+1 : 

cos0 w+1 =cosc*y n cos w +sinay n sin0 n cos0r&). 



(5.21) 



(5.22) 



The geometrical derivation goes as follows. Let Q n+l be the polar deflection and Xn+i the 
associated azimuthal deflection in the laboratory system, resulting from the n+lst collision. 
xw+ i is assumed to be measured with respect to a plane that contains the £-axis and the direction 
of motion of the neutron before the w-flst collision.) The various deflection angles and 
obliquities form a spherical triangle that is sketched in figure 5. The application of the law 
of cosines to this triangle yields the relation 



cos0 n+ i = cos0 n+1 cos0 n +sin0 w+1 sin0 n cosx»+i. 



(5.23) 



z -DIRECTION 




NEUTRON DIRECTION 
AFTER n'th COLLISION 



NEUTRON. DIRECTION 
AFTER n+l'st COLLISION 



POINT OF lul'st 
COLLISION 

Figure 5. Spherical triangle illustrating the kinematics of a neutron scattering. 
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Now, according to (4.5) and (4.10), 6 n+1 is identical with a jn . Furthermore, it can be assumed 
that the azimuthal deflection angle Xn+i is distributed randomly between and 2tt, so that it is 
legitimate to replace eosx«+i in (5.23) by cosirfS, where £<5 is a random number. Thus (5.23) 
is in fact identical with (5.22). 

5.3. Albedo Estimates 

The estimate of the number albedo A N involves merely the determination of the fraction 
N/N of the number of histories that end in reflection. 12 To obtain an expression for the 
statistical error of such an estimate we shall be somewhat formal and assign to each sampled 
history a "score" R, which has the value 1 in case of reflection and is otherwise. Thus, R is 
a random variable which assumes the value 1 with probability A N and the value with proba- 
bility 1— A N . 

The mean value of R is 

R=A N -l + (1-A N )-0=A N (5.24) 

and its variance (defined as R 2 —R ) is 

var R=A N {l—A N ). (5.25) 

Suppose now that we sample and examine 2V histories with scores RijR 2 , . . . ,AW An esti- 
mate of the number albedo is given by the average score 



-&*+*•+ • • ■+**> =w 



A&= w (R l + R2+ ■ ■ - + />V ) =xr (5-26) 



The variance of this estimate 13 is 



YarA$=jyV*rR=jj-A„(l—A N ). (5.27) 

An estimate of the energy albedo can be obtained by multiplying the number albedo 
estimate N/N hy the quantity E/E , where E is the average energy of the reflected neutrons. 
More formally, we assign to each history a score T which is equal to E/E in case of reflection 
and which vanishes otherwise. The random variable T has the mean value 

T=A N ^ (5.28) 

and a variance 

var T=T'- T*— (A N E*-A 2 N E 2 ). (5.29) 

An estimate of the energy albedo is provided by the average score 

,lI=^(T 1 + T 2+ ... + r„ o )=||. (5.30) 

The variance of this estimate is 

var Ai~j£^ {A N W-A%E 2 ) 

= WW i ^ var R+Am var E) - (5 - 31) 



12 Crude estimates of the differential albedo can be obtained by counting the neutrons that are reflected with energies and obliquities in speci- 
fied energy-angle intervals. The statistical accuracy of such estimates is given by formulas analogous to those for the number albedo. 

u Actually one does not know the exact value of An, but an estimate of var (^4^) is obtained through replacing A n in (5.27) by A N , and No 
by No— 1 in order to make the estimate unbiased. 
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The two terms contributing to this variance are both non-negative. The first arises from 
fluctuations of the number of reflected neutrons and the second from variations of the energies 
of the reflected neutrons. 

The significance of the variances var A N and var A E must be interpreted with reference 
to the probability distributions of A N and A E . Both of these quantities are the averages of 
iVo identically distributed independent random variables. According to the central limit 
theorem of statistics [14] their distributions approach a Gaussian distribution in the limit of 
large iVo- 14 One can therefore make the following statements about the deviations of the 
albedo estimates from the true albedo values : 

Probability { | A$—A N \ >kJva.TA$} = Probability { \Ag—A B \ >fcV™rZ|} =J- i e~ x2,2 dx. 

(5.32) 
5.4. Correlated Random Sampling 

To establish the dependence of the albedo on the assumed conditions (source energy and 
obliquity, collision cross sections, etc.), it is better to use a difference method than to make 
two or more separate calculations. The difference method can be designed to eliminate 
irrelevant statistical fluctuations, by keeping constant all random elements in the calculation 
which are not affected by the change of assumed conditions. Thereby greater accuracy can 
be achieved with a given amount of computation. 

We recall from section 5.2. that the following random numbers specify the nature and 
outcome of a collision: £ a determines the quantity s n +i/\(E n ) which is the distance (in mean 
free paths) traveled by the neutron between the previous collision and the collision under 
consideration; £3 determines the type of collision; £ 7 , and, for anisotropic elastic scattering 
from oxygen, additional numbers £ 7 n , £ r 2) , . . . determines the center-of -mass deflection Q c ,n+i m , 
finally, & determines the azimuthal scattering deflection x«+i- The above characteristics, 
together with the initial conditions, specify a history completely. 

Suppose now that we want to make a direct determination of the difference in the albedo 
brought about by a change of the obliquity of the incident radiation. We can do this by 
generating simultaneously two correlated sets of histories, one set with one initial obliquity, 
the second set with another. The sets should be correlated in the sense that for each history 
in one set there exists a counterpart in the second set in which exactly the same random numbers 
are used for determining the outcomes of corresponding collisions. Thus two matched 
histories are alike, except that one of them is "rotated with respect to the other' , and may 
be shorter because of the different initial obliquity. We call this case one of maximal cor- 
relation. 15 

In order to find out how a change in the assumed cross sections affects the albedo, we 
may again generate two (or more) correlated sets of histories, using identical random numbers 
for corresponding choices. Then two matched histories will at least be identical in regard to 
the distances (in mean free paths) between successive collisions, and in regard to successive 
azimuthal deflections. However, energy losses in corresponding collisions may be different. 
Consequently all the parameters of a history which depend on the energy may also be different. 
We call this case one of partial correlation. Even in this case, the use of identical random 
numbers tends to minimize accidental differences. 

It is most economical, but not absolutely necessary, to generate several correlated sets 
of histories simultaneously. In the present work we wanted to compare the albedo under 
many different conditions, and used the neutron histories also for the solution of other prob- 
lems, so that simultaneous computation of all correlated histories was not practicable. We 



u In fact, the distribution of A* N is the No-fold convolution of a binomial distribution with itself, and is very closely approximated by a Gaussian 
even for No as low as 10. The numerical results for the energy albedo indicate that the distribution of A* E is a reasonable approximation to a 
Gaussian for iVo«100. 

» This case also arises in the comparison of the reflection of neutrons after 1, 2, . . . n collisions. 
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therefore bad to devise a technique for doing separate but correlated calculations. This was 
facilitated by the use of pseudorandom number sequences, which according to (5.3) and (5.4), 
are completely determined by a "starting value" c and a "multiplier" 5 k , and can be conven- 
iently regenerated provided we remember these two parameters. 

The number of differenl random numbers required to generate two histories is usually 
not the same, even when they are correlated. In the first place, even in the case of maximal 
correlation, the two histories will tend to differ in length. Secondly, in the case of partial 
correlation there is a statistical variability in the number of random numbers £^ 1} , ^ 2) , . . . 
that may be required to select a scattering deflection by the rejection technique, according 
to (5.17). Thus it could happen thai in spite of the use of the same pseudorandom number 
sequence (with the same c and k), different portions of this sequence would be used for 
supposedly correlated histories. This difficulty was circumvented by the use of many 
sequences of pseudorandom numbers instead of just one. A primary sequence was introduced 
with an arbitrary odd starting value c and a multiplier 5 fcl . Successive members of the primary 
sequence were then used as the starting values of secondary sequences generated with a 
multiplier 5 A \ Successive secondary sequences were used for the calculation of successive 
histories. When the need arose for random sampling by the rejection technique, t ho requisite 
random numbers £$, 1} , £^ 2) , . . . were obtained from a tertiary sequence. The tertiary sequence 
used a multiplier 5* 3 , and a starting value equal to the last member of the secondary sequence 
currently being used. 10 

We now discuss t he effect of correlated sampling on Monte ( )arlo est imates of A N —A N and 
A E —A E , where A N and .l/. ; are albedos under condition 1, and A N and A,, : the albedos under 
condition 2. Proceeding as in section 5.3, we assign to histories generated under condition 1 
the scores R and 1\ and to the histories generated under condition 2 t he scores R and T. Albedo 
differences can then be estimated as follows: 

A%-A%=j^{(R l -R l ) + {R 2 -K 2 )+ . . . + (R Nq -R Nq )}, (5.33) 

At-AS=±{(T l -f 1 ) + (T 2 -fj + . . . +(T Nq -T Nq )}. (5.34) 

The variances of these difference estimates are 

var (A%— A$)= w {var R-2p(R,R)^\SLr R var R +var /7} , ( 5.35) 



~N { 



The quantity 



var (4*-A%)=^r {var T-2p E (T, T)Vvar Tvar T+var T}. (5.36) 



/ r> r»\ RR — R'R /» ~p- s 

p(R 7 R)==—==== (5.37) 

V var R var R 



is the correlation coefficient of R and K. Similarly, 



u We have set A* i = 1 J , fc 2 =15, and fa=13. 



TT—T-T 

PE (T,T)^-^==^^= (5.38) 

V var 7 var T 
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is the correlation coefficient of T and T. For uncorrelated pairs of histories, the correlation 
coefficients vanish. The objective of correlated sampling is to reduce var (A£—A$) and var 
{A%— AS) by making the correlation coefficients positive and as large as possible. 

The magnitude of the actual correlation coefficients must, in general, be estimated from 
the sampled data, and some typical values will be given in section 7. It is possible, however, 
to put a theoretical upper limit on p(R,R). Using (5.25) and (5.37) we find that 

p (R, R) = __ R%-A„A„ => (53Q) 

tJA n (1-A n )A n (1-A n ) 

Without loss of generality we may assume that A N <A N . 
Therefore 

probability {i?=l and #=l}<probabiHty{5=l}=4^ (5.40) 

so that 

RR<A N . (5.41) 

Accordingly the correlation coefficient must obey the inequality 



The greatest possible value of p(R,R) corresponds to the case of maximal correlation discussed 
above. It leads to a reduction of var (A N —A N ) by a factor 



1 



A N | 1-^ (5.43) 

A N 1 A N 

6. Computations and Results 
6.1. Machine Computations 

The Monte Carlo calculations were programed for the IBM-704 computer at the National 
Bureau of Standards. 17 The neutron histories, generated in order to determine the albedo, 
were concurrently used for the calculation of the neutron flux in infinite media and in slabs, 18 
so that it is difficult to state how much computer memory and computing time were required 
for the albedo problem alone. Altogether, a little more than 5,000 words of magnetic core 
memory were required of which not less than one-third were allocated to the albedo problem. 
For each case treated, 3,000 neutron histories were generated and analyzed. This amount of 
computation gave the desired accuracy and required approximately 45 min of computer time. 
Not more than half of this time was required by the albedo problem. The use of magnetic tapes 
was not necessary, but was an optional device for storing the results of one machine run and 
combining them with the results of a later run. Intermediate print-outs of the results were 
made after groups of 1,000 histories had been processed, in order to follow the course of the 
computations and to have a check on the performance of the computer. 

In line with the exploratory nature of the investigation, the machine program was kept 
flexible. This made possible the exploration of new aspects of the problem as they emerged in 



17 The program was based in part on an earlier code [15] developed at the Westinghouse Atomic Power Division. The authors are indebted to 
Dr. B. II. Mount and Mrs. Heidi Kuehn who very graciously made their code available to us. 

18 To be reported in a later publication. 
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the course of the work. In the simplest version of the program only elastic scattering from 
hydrogen and oxygen was taken into account, and the print-out of results was limited to tables 
of the differentia] albedo. In the latest version all relevant collision cross sections were in- 
cluded, and much additional information was printed out that had been found necessary for an 
understanding of the physical and statistical significance of the Monte Carlo results. This 
additional information included the mean and mean square energy of the reflected neutrons, 
the contributions of different orders of scattering to the number albedo and energy albedo, 
listings of the particular histories in which reflection occurred and of the energies of the reflected 
neutrons, and a statistical analysis of tin 4 pseudorandom numbers used. Intermediate versions 
of the program yielded only some of this information. 

The analytical formulas for the single-scattering albedo, developed in section 4, were used 
mainly for the interpolation of results, and single-scattering computations were done with the 
Monte Carlo code. This was simpler to do once the machine program was available but slightly 
less accurate. 

Altogether 38 different cases were treated which are listed in table 6.1. This table contains 
the source energy and obliquity as well as a brief indication of the assumptions about the cross 
sections. The computations fall into two groups: 30 cases were done with one set of pseudo- 
random numbers (analyzed in app. B), and the remaining 8 cases with another set. A high 
degree of correlation was introduced by the repeated use of the same set of pseudorandom num- 
bers in order to increase the accuracy of intercomparisons. We believed this to be a sound 
procedure, because the set was large enough (consisting of approximately 237,000 numbers) 
to be statistically representative. 

Table 6.1. Summary of calculations 

The entry in the left column serves as index to the differential-albedo tables in appendix C. Each result for the number albedo An and the 
energy alhedo Ar is based on the analysis of 3,000 neutron histories, with SAx and Mi indieat lug the Stal Istical standard deviations. All runs 
were done with the same set of pseudorandom numbers (described in app. B), except the runs indicated hy an asterisk (*) which were all done with 
another set of pseudorandom numbers. 



Table 



1 

None 

2a 

2b 

2c 

3a 

3b 
4a 
4b 
5a 
5b 

6a 
fib 
7a 
7b 
7c 
None 



8b 
8c 
9a 
9b 

10a 
10b 
11a 
lib 
12a 
12b 

13 
14a 
14b 
15a 

15b 
15c 



16a 
16b 
16c 
16d 



A'o(Mev) 



0.3 

1 
1 
1 
1 



14 
14 
14 
14 



0° 
0° 
0° 
0° 
0° 

45° 
45° 
60° 
00° 
75° 
75° 

90° 
90° 
0° 
0° 
0° 
0° 

0° 

0° 

0° 

45° 

45° 

60° 
60° 

75° 
75° 
90° 
90° 



(t) 



A„ 



0. 296 
.370 
.361 
.381 
. 113 

. 450 
. 114 
.518 
. 142 
.619 
.213 

. 767 
.493 
.061 
. 131 
.163 
.168 

.178 
.167 
.028 
.232 
.035 

.315 
. 066 
.472 
.143 

.728 
.480 

. 371 
.157 
. 171 
. 107 

.138 
.131 



.124 
.163 
.165 
.064 



M« 



0. 008 
.000 

.009 
.009 
.006 

. 009 

.000 
.009 
.006 
.009 
. 007 

.008 
.009 
.004 
.006 
.007 
.007 

.007 
.007 
.003 
.008 
.003 

. 008 
.005 
.009 
. 006 
. 008 
.009 

. 1)09 
.007 
. 007 
.006 

.006 
.006 



.006 
.007 
.007 
.004 



0.048 
(. 135) 
(.115) 
.148 
. 091 

. 159 
.085 
. 183 
.099 
.241 
.144 

.415 
.347 
.023 
.038 
.044 
(. 044) 

(. 049) 
.041 
.023 
.052 
.023 

(. 080) 
(. 039) 
.144 
.084 
(. 339) 
(. 290) 

(.110) 
.052 
.059 
. 027 

.041 
.047 



.040 
.038 
.079 
.023 



0. 005 



. 003 
. 004 
.002 

.003 



.003 
.002 
.004 
.002 



Assumptions 



(*). 
(*). 

(*); elastic scattering from oxygen assumed isotropic. 

Single scattering. 

Single scattering. 
Single scattering. 
Single scattering. 

Single scattering. 

2-cm slab; see section 7.9; (*). 

5-cm slab; see section 7.9; (*). 

10-cin Slab; see section 7.9; (*). 

(*). 

(*); elastic scattering from oxygen assumed isotropic. 

Single scattering. 

Single scattering. 

Single scattering. 
Single scattering. 
Single scattering. 



Absorption disregarded. 

Excitation of 6.1-Mev level in oxygen by inelastic scatter- 
ing. 

Absorption disregarded. 

A bsorption disregarded. Inelastic scattering from oxygen 
treated as elastic. 

Excitation of 6.1-Mcv level in oxygen by inelastic scatter- 
ing. 

Elastic scattering from oxygen assumed isotropic; excita- 
tion of 6.1-Mev level in oxygen by inelastic scattering. 

Absorption disregarded. Inelastic scattering from oxygen 
treated as elastic. 

Inelastic scattering treated as absorption. 



(t) Isotropic source. 
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6.2. Differential Albedo 

Tables 1 to 16d of appendix C contain the Monte Carlo results for the differential albedo 
in ( he lornj in which they were printed out by the computer. The number of reflected neutrons 
(for N — 3,000 neutrons) is given in a two-way classification according to the spectral energy, 
E, and the cosine of the obliquity of emergence, cos 6. By normalizing the entries in the tables, 
i.e., by dividing through by N , one obtains an estimate of the differential albedo averaged over 
a small energy-angle interval. Unnormalized numbers are given in the tables in order to pre- 
serve their statistical significance, the standard deviation of each eiitr}^ being proportional to its 
square root. The energy classification is divided into 10 intervals, labeled i=l } 2, . . . 10, such 
that the iih interval pertains to neutrons with energies between (ll—i)/\0E o and (10— i)/10E . 
The last interval extends from E /10 to 0.5 ev rather than to zero energy. The obliquity 
classification is also divided into 10 intervals, labeled ifc=l, 2, . . . , 10, such that the kth 
interval pertains to neutrons with obliquity cosines between (1— k)/10 and —k/10. Two sets 
of marginal totals are also given which represent the angular distribution integrated over all 
energies, and the spectral distribution integrated overall directions, respectively. 

6.3. Number Albedo and Energy Albedo 

Table 6.1 contains all the results obtained for the number albedo A N and the energy albedo 
A E . 19 The number albedo was obtained by adding the marginal totals in the corresponding 
differential-albedo tables and dividing the result bj N = 3,000. The energy albedo was obtained 
as the product of the number albedo and the average energy E of the reflected neutrons. In 
most of the calculations A E was computed directly in the course of the machine calculations. 
In some of the earlier cases, E had to be estimated from the differential-albedo tables. 20 Nu- 
merical experimentation indicated that such estimates could be made with an accuracy of 
approximately 1 percent. Values of the energy albedo based on such estimates of E are indicated 
in table 6.1. by being placed in parentheses. 

The standard deviations 5A n =-\v2lt A n and bA E =-\v&Y A E are also shown in table 6.1- 
They were computed according to (5.27) and (5.31), respectively. The determination of 5A E 
required knowledge of E as well as of E 2 . Accordingly, 8A E is given only for those cases in 
which E 2 was obtained as direct output of the machine program. The relative standard devia- 
tion 5A N /A N , obtained with 3,000 histories, is of the order of 0.05. The corresponding relative 
standard deviation 5A E /A E , in the cases examined and presumably also in the other cases, is 
roughly \){ times greater than bA N jA N . 

7. Analysis of the Results 
7.1. Single-Scattering Contribution to the Albedo 

Table 7.1 contains a summary of theratios^4 ( ^ ) /-4ivand-4 ( ^ ) /^f or various cases. In general, 
single scattering contributes about 20 percent of the number albedo and more than 50 percent 
of the energy albedo. The contributions are particularly large at 1 Mev, because at this energy 
the ratio of the scattering cross section of oxygen to that of hydrogen is large, and the differential 
ox3^gen scattering cross section has a peak in the backward direction. 



1& In sec. 5.3. the albedo estimates were indicated by A*s and A*b in order to distinguish them from the true values An and Ae. In the 
following, we shall ignore this distinction and simply use An and As to label the Monte Carlo results. 

2 o The most accurate estimation procedure was found to be one based on the "cumulative albedo , ' I dE' I sindd9A(E' ,d;Eo,do) , which can 



be obtained from the differential-albedo tables for Ei=i/10E , i=l, 2. . . . 10, and by interpolation for other values of E. It was also necessary to 
take into account the characteristics of the spectrum of the single-scattering albedo, in order to obtain good accuracy. 



120 



The estimation of the albedo ratios is a problem with maximal correlation as defined in 
section o.4. We have estimated the standard deviations of these ratios in a few typical cases 
and find. For example, that A { ^/A N =0M±0.03 and , V,? /A E =0M±0.02 for E =l Mev 
and O =O°. 

The ratio A E } /A B is a relatively slowly varying function of the source energy, whereas A E l) 
and A E separately vary rapidly. This suggests that the energy albedo for arbitrary energies 
eonld be obtained with fair accuracy by combining a single-scattering calculation with an 
interpolated value of A k E \A E . 

Table 7.1. Ratio of single-scattering albedo to total albedo. (The results for 6, 9, and 14 Mev correspond to cases 

l^a, 15a, and 16a, respectively, as listed in tubh H I 
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7.2. Contribution of Higher Orders ol Scattering 

Figure 6 shows die contributions of successive orders of scattering to the albedo for a few 
typical cases. As many as twenty orders are shown to contribute to the number albedo, 21 
whereas only about five orders contribute significantly to the energy albedo. 'Flic relative 
importance of the first few orders of scattering increases with the source obliquity and energy. 
In general, an analytical orders-of-scattering approach does not seem to provide a practical 
alternative to a Monte Carlo calculation. The contribution of the second order of scattering 
is not sufficient for a precise albedo determination, and calculations for the higher orders would 
be prohibitively complex. 

7.3. Energy Spectra of Reflected Neutrons 

Figure 7 contains histograms of the energy distribution of reflected neutrons for incident 
1-Mev beams with various obliquities. The contributions of single scattering and multiple 
scattering to the spectra are indicated. The spectra are characterized by two rather sharp 
peaks: one, primarily due to single scattering, occurs at the upper end of the spectra ; tbe other, 
due to multiple scattering, at the low-energy end. The single-scattering contribution can in 
turn be resolved into two components, due to scattering from oxygen and hydrogen. For 
perpendicular incidence, only single-scattering reflection from oxygen occurs. As the source 
obliquity is increased, hydrogen scattering becomes possible, and the oxygen contribution is 
shifted toward higher energies. At 45° and 60° the oxygen and hydrogen components are 
clearly separated, whereas at 75° and 90° they merge. In the limiting case of grazing incidence 
there is a hydrogen component distributed uniformly over all energies, and an oxygen component 
which gives rise to a peak at an energy slightly below the source energy. Further details about 
the oxygen single-scattering contribution can be obtained from appendix B (fig. 15). The 
energy spectra for other cases are similar to those shown in figure 7 and can be obtained from 
the tables in appendix C. 



21 Some contributions from orders beyond the twentieth would have occurred if our calculations had not been cut off at an energy of 0.5 ev. 
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Figure 6. Contributions of successive orders of scattering to the albedo. 
7.4. Angular Distribution of Reflected Energy 
Figure 8 shows plots, in histogram form, of the quantity 

A E (d; Eo, Bo) = A(E, 6; E , 6 )§dE, 
Jo A) 



(7.1) 



for E =S Mev and various values of O . 22 The angular distribution is dominated by the single- 
scattering contribution. In particular, it should be noted that for large source obliquities the 
distribution does not follow the cosine-law which holds for perpendicular incidence and has 
often been assumed to be universally valid. Rather the emerging radiation tends to peak for 
cos0=O; i.e., in directions parallel to the boundary of the semi-infinite medium. This is due to 

diftJ^H^ °y i l e an g ular distribution, because each angular element was obtained by multiplying the 

mtervals elemeDt of the corresponding table m app. C by the midpoint of the appropriate energy interval and then summing over-all energy 
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1.0 



Figure 7. Energy spectra of reflected neutrons. 
The shaded areas indicate the contribution of single-scattering. Source energy 1 Mev. 

the behavior of the single-scattering component. For Z? =3 Mev and O =9O°, we find from 
eq (A.21), (A. 22), and (A. 24) of appendix A, and from the cross sections given in figure 2, that 

^^-0.049(1.00+1.05 sin 2 0+7.14 sin 3 0), (7.2) 

which is in agreement with the Monte Carlo calculations. 

7.5. Dependence of Albedo on Source Energy and Obliquity 

The dependence on the source energy is irregular. It is correlated with the complicated 
energy-dependence of the attenuation coefficient for water, and is further influenced by the 
relative probabilities for scattering by hydrogen and oxygen, and by the energy-dependent 
anisotropy of elastic scattering from oxygen. Just like tin 4 attenuation coefficient, the number 
albedo and energy albedo have a peak at 1 Mev where they are 2 to 3 times larger than at the 
other energies considered. The correlation between attenuation coefficients and albedo is 
complicated, however, by the interaction of the various factors, so that interpolation with 
respect to source energy is not easy, and is perhaps done best by the single-scattering analysis 
indicated in section 7.1. 23 

23 For the analysis of an albedo experiment with a broad source spectrum (such as the Po-Be source used by Strickler, Gilbert, and Auxier 
[16]), it might be preferable to make a Monte Carlo calculation with a realistic source spectrum instead of combining the results of calculations for 
monoenergetic sources. 
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Figure 8. Angular distribution of reflected energy. 
Source energy 3 Mev. A single-scattering curve for 0o=9O° is shown which corresponds to eq (7.2). 

The dependence on the source obliquity is smooth and interpolation should be easy. By 
way of example, figure 9 contains curves of A N and A E as functions of d , for source energies of 1 
and 3 Mev. All four curves shown have similar shapes, starting out with a rather flat slope 
near O =O° and rising very steeply near 6 =90°. The closer 6 is to grazing incidence, the 
smaller are the differences of the albedos for different source energies. 

7.6. Effect of Absorption and Inelastic Scattering 

If absorption were disregarded, the number albedo at 6 Mev would be increased b}^ 9 per- 
cent and the energy albedo by 13 percent. At 9 Mev, the corresponding increases would be 29 
percent and 52 percent, respectively. The effect of absorption is greater on the energy albedo 
than on the number albedo because the relative probability of absorption is an increasing func- 
tion of the energy (up to 14 Mev). The effect of absorption on the energy spectrum is shown in 
figure 10. Although the effect is large, it introduces no large uncertainty into the calculations 
because the cross section for n — a and n—p processes are reasonably well-known. 

The effect of inelastic scattering is more difficult to determine because so little is known 
about the cross section for this process. We have made computations under three different 
assumptions: (1) Absorption taken into account, inelastic scattering assumed to excite oxygen 
to a level 6.1 Mev above the ground state; (2) absorption taken into account, inelastic scattering 
also treated as absorption; (3) absorption disregarded, inelastic scattering treated as if it were 
elastic. 

As has been discussed in section 3.2, assumption (1) is presumably the most realistic. 
Assumption (3) is the least realistic and has been included mainly because it has often been made 
in previous neutron transport calculations. Assumptions (1) and (2) are limiting assumptions 
in the sense that they allow for the minimum and maximum possible energy loss that a neutron 
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Figure 10. Reduction of energy spectrum due to absorption, indicated by shaded areas. 
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could suffer in an inelastic collision. The albedo calculated on assumption (3) is almost twice 
as great as that calculated on the realistic assumption (1), whereas the albedo calculated on 
assumption (2) is only about half as large. These extreme variations point up the need for 
accurate knowledge of the inelastic scattering cross section. 

The energy spectra obtained with a 14-Mev source under the three assumptions are shown 
in figure 11. The spectrum in case (1) shows three peaks, due to elastic single scattering, in- 
elastic single scattering, and multiple scattering. The inelastic scattering peak disappears in 
case (2) and merges with the peak due to elastic single scattering in case (3). More accurate 
treatment of inelastic scattering would probably give a spectrum similar to that of case (1), but 
with the middle peak reduced in size and "smeared out" over the lower part of the spectrum. 



A(E) 






Figure 11. Effect of inelastic scattering on the energy spectrum. 

Source energy 14 Mev. Histograms (i), (ii), and (iii) correspond to the different assumptions in sec. 7.6, and were obtained from tables 16a 
16d, and 16c of app. C, respectively. Perpendicular incidence. 



7.7. Effect of the Anisotropy of Elastic Scattering 

The differential elastic scattering cross section is not well-known for most elements. It is 
therefore of interest to know how much the albedo would be changed by a change of the assumed 
distribution of elastic scattering deflections. To explore this situation, we have compared the 
albedo calculated with the best available cross sections for oxygen (reference [9], see also fig. 3) 
with the albedo calculated on the assumption that the elastic scattering deflections from oxygen 
are isotropic in the center-of-mass system. The difference between the two assumptions is 
greater than any expected realistic modification of the oxygen cross sections. Comparisons 
were made for source energies at 1, 3, and 14 Mev. At 1 Mev, A N is decreased by 2.5 percent 
and A E by 11 percent. Finally, at 14 Mev, A N is increased by 31 percent, whereas A E is de- 
creased by 5 percent. In figure 12 the spectra obtained under the two different assumptions 
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Figure 12. Effect of anisotropy of elastic scattering from oxygen on the energy spectrum. 

Perpendicular incidence. 



are compared. There is relatively little change, the main effect being a shift of the single- 
scattering peak toward higher energies in the case of isotropic scattering. 

The albedo differences quoted above were obtained by correlated sampling (case of partial 
correlation). A listing of all significant statistical parameters for the comparison at 14 Mev is 
given in table 7.2 from which it appears that by means of correlated sampling an accuracy was 
achieved that would have required a sample size 3 to 4 times as large with ordinary sampling. 



Table 7.2. Estimate of albedo difference by correlated sampling 
Source energy 14 Mev, perpendicular incidence. Elastic scattering from oxygen assumed to be anisotropic in case 16a, isotropic in case 16b. 



Number albedo: 
An, case 16a, 
An, case 16b, 
<4#— Jjv=0.039=fcO.0O5 (with correlated sampling), 

±0.009 (without correlated sampling), 
Correlated sampling increases effective sample size by factor 2.9, 
Correlation coefficient p =0.664, 
Maximal correlation coefficient p m ax =0.853, 



Energy albedo: 
As, case 16a, 
Ae, case 16b, 
vljr-^4if=-0.0012±0.0016 (with correlated sampling), 

±0.0033 (without correlated sampling), 
Correlated sampling increases effective sample size by factor 3.1 
Correlation coefficient pj?=0.771. 
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7.8. Low-Energy End of the Spectrum 

Details of the low-energy part of the spectrum were obtained as a bj^product of the detailed 
information obtained for the statistical analysis of the correlated calculations at 14 Mev. 
Figure 13 shows a plot of the spectrum on a logarithmic scale which emphasizes the low energies. 
The spectral shape at low energies is quite similar to that obtained by Spinney [5] who calculated 
by age theory the albedo for 2-Mev neutrons incident on concrete. This suggests that com- 
putations of the spectrum down to very low energies could advantageously be done by a 
combination of Monte Carlo and age theory, with Monte Carlo covering the energy region 
not very far below the source energy in which age theor}^ does not apply. 
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Figure 13. Energy spectrum of reflected neutrons on semilog scale. 
Source energy 14 Mev, perpendicular incidence. 



7.9. Slab Albedo 

The deeper a neutron penetrates into the medium, the smaller is its chance of returning. 
Therefore, the collisions leading to reflection take place rather close to the surface of the 
medium, and thin slabs have nearly the same reflecting power as a semi-infinite medium. 
This is shown by the results in figure 14 pertaining to the reflection of 3-Mev neutrons incident 
on slabs with thicknesses of 2, 5, 10, and 20 cm. (These thicknesses are to be compared with 
the transport mean-free path in water at 3 Mev which is 5.1 cm.) The dependence of the 
number albedo on the slab thickness is similar to that found by Foderaro and Obenshain [6] 
at 0.89 Mev. The energy albedo reaches a limiting value at an even smaller slab thickness 
than the number albedo, because the energy degradation in successive collisions makes it 
necessary for a neutron to be reflected very early in its history if it is to carry back an appre- 
ciable amount of energy. 

128 




Figure 14. Dependence of slab albedo on dub thickness. 
Source energy :* Mev, perpendicular incidence, a. Number and energy albedo and b. Energy spectrum. 

| | 2-cm slab 
^-f-Q 5-cm slab 
^+Ey + Q semi-infinite medium 

8. Appendix A. Integrals Over the Differential Single-Scattering Albedo 

8.1. Energy Spectrum 

The evaluation of (4.15) leads to the following expression for the energy spectrum of neu- 
trons reflected from a semi-infinite medium after exactly one collision: 



where 



A™ (E;EM = i:vj(Eo)F j (E;E ) { U{a h b 3 ) -kV(aj,bj) } , 
;=i 



a^-=cos (aj-\-6o) 
bj=eos (a j — O ) 
Jc=(\q/\) cos d , 
and where U(dj,bj) and V(a j} bj) are defined as follows: 

(1) If E [ VM|^-co^g +sin g .J< K <^ 

U(a J ,b J ) = V(a j ,b J )=0 

(2) If Eo |" VM^-co^o- S in e - J< E < E(> ^/M 



jKj— cos 2 o +sin d 



Mj+1 



J 



Uia^-lsin-^ 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
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iik<Cbj, V(a j} b J )=— F 



^(bj—k)(k—aj) 



{log (bj— a } )k— log [{b 1 -\-a j )k—2a j bj 



—2^jajbj(bj—k) (a,—k)]} ; 



.„ x , w u 1 fl 1 • _ i r(b J +a J )k— 2ajb~\\ 

iik>Oj, F(a y ,6 < )= , J „ sin ' — ., — ., >-• 

" ' ^(k-bj)(k-aj)\2 * L (bj-a^k J J 



^ (A.5) 



£7(a,A-) = l 



F(a y ,6 y ) = 



1 



(4) If E<E 



L M,+l J 



■\l(k—bj)(k—aj) 



U(a J ,b J ) = V(a ] ,b j )=0. 



(A.6) 
(A .7) 

(A.8) 



In certain limiting cases, the formula for the energy spectrum becomes much simpler than 
the above expressions. For perpendicular incidence, the energy must lie in the range 



T M,g,-l l 



and the spectrum is 



i 



A«KE; E 0j O)=j:p j (Eo)F j (E; E ) 



j=i 



1 



X cos otj 



For grazing incidence (0 O — ^tt/2), the energy range is 



E { 



and the spectrum is 



, r M,g,-i 7 F< r M,g,+i T 



^<" (#; So, |)=gpX£o) J Fj (E; E c 



(A.9) 
(A.10) 

(A.ll) 
(A.12) 



Figure 15 shows energy spectra of neutrons reflected after a single scattering from a semi- 
infinite medium of pure oxygen (E =l Mev, 00=0°, 45°, and 90°). These spectra were 
obtained by numerical evaluation of the factor 

F 2 (E;E ){U(a 2 ,b 2 )-kV(a 2 ,b 2 )} 

in (A.l). They are characterized by a peak at the low-energy limit which is rather independent 
of the source obliquity 6 . In addition there is another peak at the high-energy limit in the 
case of grazing incidence. 

For perpendicular and grazing incidence the spectrum for reflection from water can be 
derived very simply from the corresponding spectrum for oxygen. For 6 =0°, reflection from 
hydrogen is impossible after one scattering, and it is sufficient to renormalize the oxgyen 
spectrum in figure 15, multiplying it by p 2 (E ) =0.485. For O =9O O , one must not only normal- 
ize the oxygen contribution but also add a hydrogen contribution which is energy-independent 
and has the value p 1 (E ) 1 AF 1 (E' } E ) = 0.515 (1/2E ) = 0.2575. 
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Figure 15. Energy spectrum of neutrons reflected after one collision. 
Source energy 1 Mev, oxygen medium. 



8.2. Angular Distribution 



The angular distribution of neutrons reflected after one collision is even more complicated 
than their energy spectrum. We shall confine ourselves to the discussion of some limiting 
cases. The characteristics of the angular distribution are expressible in terms of the integrals 



J cos0— (X /X) cos#o 



(A.13) 



which are needed for the evaluation of (4.16). The angular distribution of reflected energy, 
A™ (9;E Q ,d ), depends on the integrals 



J cos 6 — (Xo/X) cos O 



(A.14) 



The integrations in (A.13) and (A.14) are facilitated by changing the variable of integration 
from E to cos a,, noting that 

„_ w f eos aj+Vcos 2 aj+M jK'j- 1"| 2 
E - E ±- ~^+i J 

and that 

dE 2E Q f j lKJ2 T 2 ^ — cos 2 <x, ^i 

1 = /*>r i ix2 -< 2 cos aj+^MjKj— sm 2 aj+-j= j =_ V- 



(A.15) 
(A.16) 



We now take up limiting cases, and drop the index j for simplicity. 

Case (1): Heavy nucleus, M^> 1 

Perpendicular incidence, O = O° 
Elastic scattering, K= 1 . 
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Note that 

dE J1E 
d cos a M J Y (A.17) 

c «a, F(E;E )~~f(a;E ) 
Substituting expression (4.13) for H(E,6; E ,0), we find that 

K{e,E Q fl)~ co ^__ 1 , (A.18) 

and that 

h E (d;E ,0) ~E h(d;E ,0) (A.19) 

Case (2): Heavy nucleus, M^>1 

Grazing incidence , 6 =90° 
Elastic scattering, K=l. 

The approximations (A.17) may again be used. For H( E,d;E -J we substitute expression 
(4.14), the result being 

\ z / ^ J-sia&ysm^— cos 2 a 

Assuming for J(a;E ) a Legendre expansion of the type (3.9), with expansion coefficients 
jl (E ), we can integrate (A.20) and find that 

fc(0;#o,f)=|j^/ 2 /(So) (-l) e ^~,^(cos 0), (A.21) 

(8;Eo,l)~Eoh(e;E ,iy (A.22) 



and that 



Case (3): Hydrogen nucleus, M=l 
Grazing incidence, # =90°. 



We have j(e c ;E ) = 1/2 and 0<c*<7r/2. Using (A.14), (A.15), and (A.16), we find that 

i(*w T \ 2f sin « cosadcosa 2 . /Aoo , 

M 0;#o,s )=- ,^ 2/l 2 - = -sm^, (A.23) 

\ ^/ ttJo Vsiir0— cos 2 a T 



Vsin 2 0— cos 2 a: T 

and that 

, / fl F ir\ 2£ , f sin * cos 3 adcosa 4E . 3 , , A 0yI . 

M 0;£o,^ )=— - , =-^- sm 3 0. (A.24) 

\ -/ * Jo ysm 2 (9— cos 2 a * 

9. Appendix B. Pseudorandom Numbers Generated by Congruential 

Multiplication 

9. 1 . Proof of Their Uniform Distribution 

We consider the pseudorandom number sequence {? w =2" 35 x w } (n=0,l, . . ., 2 33 — 1) 
derived from a sequence of integers {x n } which is generated according to (5.4) or by the 
equivalent rule 

x n =c5 kn mod 2 35 , (B.l) 
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where c and k are odd integers. The sequence {£„} consists of 2 33 different numbers which are 
distributed with uniform density over the interval (0,1). This can be proved with the use of 
standard Dumber-theoretical results [17]. However, Dr. Morris Newman 24 has given a simple 
proof proceeding from first principles which we reproduce here in a somewhat amplified form. 
The uniformity of the distribution of the pseudorandom numbers obtained by congruential 
multiplication is a consequence of the following theorem. 

Theorem: The sequence of integers 

{x , x h x 2 , . . ., a^-l} (B.2) 

<k nerated according to (B.l), is a permutation oj the sequence 

{1, 5, 9, . . ., 2 35 -3} (B.3) 

or of the sequence 

{3, 7, 11, . . ., 2 35 -l} (B.4) 

depending on whether the starting number x Q =c has remainder 1 or 3 when divided by ./. 
The proof of this theorem is based on the following lemma : 
The integer 5 k2M —l is divisible by 2 M+2 but by no higher power of 2, if k is an odd integer. 

To see this we consider the factorization 

5«^— l = (5«*- 1 +l)(6^ jr " x — l)==(5« ir " 1 +l)(5* 2lf " a +l) . . . (5*+l)(5*— 1). (B.5) 

The factor 

5*-l = (l+4)*-l=4HSf fc V J ( B - 6 ) 

is divisible by 4 but by no higher power of 2. The other factors in (B.5) have the form 

5 M« +1 = ( 1+ 4)*2- +1=2 + W A ? W,0<m<M-l (B.7) 

;=i \J / 

and are each divisible by 2 but by no higher power of 2. By combining all the factors we obtain 
the lemma. 

We now consider the sequence {c5 kn }, n=0,l, . . . 2 33 — 1. The difference between any 
two terms in this sequence, say c5 k{n+nn and c5 kn , can be expressed as 

Atn,n')=c5 kn (5 kn '-1). (B.8) 

We may set 

n' = 2% (B.9) 

where a is an integer <C2 33 , and where b is an odd integer. By application of the lemma it 
follows that 5**'— 1, and therefore A(n,n'), is divisible only by powers of 2 smaller than 2 35 . 
This implies that the integers x n in the sequence (B.2) are all different from each other. 

It can be seen from (B.8) that the differences between the x n 's are all multiples of 4. Thus 
two different sequences {x n } can be formed, one of which is a permutation of (B.3) and the other 
of (B.4), depending on the starting number c. Each of the two sequences contains 2 33 num- 
bers, and there is no overlap between them, so that together they exhaust the 2 34 odd integers 
<2 35 . This completes the proof of the theorem. Finally, we note that by computing the 
x n 's modulo 2 N rather than modulo 2 35 , and by using a starting number c, that is not odd but 
contains a factor 2 7 , one would obtain a sequence of 2 N ~ 2 ~ 7 different numbers. 



24 Informal memorandum. Applied Mathematics Laboratory, NBS, Oct. 27, 1951. 
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9.2. Statistical Indications of Their Randomness 



In order to be acceptable as a source of random numbers for a Monte Carlo calculation, 
the pseudorandom number sequence {% n } must fulfill the following requirements: 

(1) In any subsequence of reasonable length, such as would be used in generating neutron 
histories, the £ n 's must be distributed with uniform density in the interval (0,1). 

(2) Successive £ w 's must be uncorrelated. 

These two conditions cannot be satisfied rigorously, but empirical statistical tests indicate 
that they are satisfied in very good approximation. For example, elaborate tests on 16,384 
pseudorandom numbers have been reported in the memorandum mentioned in footnote 24. 

We have made some rough-and-ready statistical tests on 237,000 pseudorandom numbers 
that were used to generate a set of 3,000 neutron histories. 25 These numbers were divided 
into three sets, corresponding to the three types of pseudorandom number sequences (primary, 
secondary, and tertiary) described in section 5.4. Set P includes the numbers in the primary 
sequence, set S the numbers in all the secondary sequences, and set T the numbers in all the 
tertiary sequences. Separate tests were applied to each set. In order to gauge the uni- 
formity of the distribution of the numbers, we computed the averages % n and £ 2 . In order to 
check on the lack of correlation between successive pseudorandom numbers we computed the 
mean square difference of successive numbers, (£ n+i — £ n ) 2 . As a further test of possible correla- 
tion, the algebraic sign of the differences between successive pseudorandom numbers was 
examined, and an index t was computed which was defined as the excess of positive over nega- 
tive signs. Table B.l compares the actual values of all these quantities with the values that 
one would expect if the pseudorandom numbers were truly random. We infer from the table 
that the pseudorandom numbers are sufficiently random for our purpose. 

Table B.l. Analysis of the randomness of the pseudorandom numbers 

In order to indicate the significance of the deviations of the empirical from the expected values of the various statistics, the standard devia- 
tions of the empirical values are given. The pseudorandom numbers were obtained by starting a primary sequence with the initial value 
c= 16,743,785,845. The number of pseudorandom numbers in each set is indicated within brackets. 



Statistic 


Expected 
value 


Empirical value 




Set P [3,000] 


Set S [220,776] 


Set T [13,384] 


IT. 


H 
H 
H 




0. 49402±0. 00527 
.32937± .00.544 
.16743± .00360 

28 ±55 


0. 49955±0. 00061 
.33359± .00064 
. 16579± . 00042 
-638 ±470 


0. 50233±0. 00250 


il 


33580± 00258 


(£ n +l— In) 2 - 


.16461± .00171 


t 


34 ±116 







10. Appendix C. Tables of Differential Albedo 

The following tables are indexed by table 6.1 and described in section 6.2 of the text. 















Table 1 


.. 












cos 




\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 


3 
1 


7 
2 


4 

1 
2 


8 
3 
1 


12 
1 
2 


14 

5 


14 
1 
1 


6 
3 
2 
3 
2 
2 
2 
6 
101 






68 




3... 


15 
3 

5 
3 
2 
4 

124 


18 
3 
3 
7 
3 
2 
6 
142 


50 




4 


14 


E 


5 






3 

5 
4 
1 
5 
75 


10 




6 










2 
1 
2 

4 

48 


4 

5 
3 
6 
91 


25 




7 . 










18 




8. . 








1 
2 

49 


13 




9... 


1 




1 
16 


35 




10 


10 


656 












Sum 


5 


19 


24 


64 


72 


112 


125 


127 


157 


184 


889 



































25 This set of numbers was used repeatedly for most of our calculations. 
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Table 2a. 



cos e 





\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 






1 

Mi 














1 




2 


5 


10 


20 
7 
2 
1 
3 
1 
1 
3 

43 


32 

4 
2 
5 
2 
5 
1 
5 
52 


37 
7 
4 
3 
4 
9 
5 
4 

48 


48 
4 
3 

5 

8 

5 

6 

52 


30 
14 
3 
4 
2 
3 

15 
94 


3 
40 
3 
2 
11 
9 
3 
11 
102 




201 




3 


41 
5 
5 
5 
6 
10 
10 
123 


117 




4 


1 


1 
2 
1 


2 
1 


26 


E 


5. 


28 




6 




33 




7— 




2 
2 


43 




8 




1 
1 
10 


35 




9 


1 
3 


56 




10 


17 


544 










Sum . 


10 


26 


41 


81 


108 


121 


136 


172 


184 


205 


1,084 









Table 2b. 





'^\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 . 

























2 


1 

1 




5 
2 


10 
5 
2 


21 
6 

3 

7 

5 

4 

49 


32 
5 
3 
5 
4 
5 

10 
5 

51 


46 
12 
3 

4 

1 

2 

11 

51 


56 

34 
6. 

11 
8 
4 
5 
5 

64 


1 
86 

8 
3 
7 
7 
8 
9 
86 




172 




3... 


2 


113 
12 

7 
7 
4 
11 
10 
88 


266 




4. _ 


39 


E 


5 




1 
1 
1 




35 




6 




2 




34 




7 




3 

8 
4 
37 


32 




8 




3 


52 




9. _ 


2 
4 


5 

7 


55 




10 


22 


459 










Sum 


8 


17 


34 


69 


101 


120 


135 


193 


215 


252 


1,144 









Table 2c. 





^\ k 

i \ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 


1 




4 


9 


19 


31 


42 


56 

21 






162 




3 




68 


88 


177 




4 



















E 


5 



























6 



























7 



























8 



























9 



























10 





















































Sum 


1 





4 


9 


19 


31 


42 


77 


68 


88 


339 









Table 3a. 





i ^X 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1... 

2 

3 

4 _ 


2 

8 


2 

4 
11 

1 
4 


1 

13 
14 

2 

1 
2 
2 
9 
9 
24 


3 

9 

18 
7 
3 
2 
3 
6 
7 
42 


2 

Hi 

32 
7 
4 
2 
5 

12 
3 

56 












9 




11 

21 
6 

10 
7 
2 

10 
9 

62 


17 
33 

7 
9 
8 
5 
7 
10 
59 


21 
49 
10 
2 
7 
3 
6 
12 
84 


18 
47 
6 
4 
10 
9 
14 
10 
100 


37 

45 
8 
4 
7 
4 
9 
19 

123 


150 
271 
55 


/•' 


5... 


42 




6 


46 




7 


10 
4 
5 

11 


44 




8 


78 




9 


86 




10 .. 


569 










Sum ._ 


21 


52 


77 


100 


139 


138 


155 


194 


218 


256 


1,350 









135 



Table 3b. 





\ k 


1 


2 


3 


4 


5 


G 


7 


8 


9 


10 


Sum 




1 

2 

3 


1 
4 


2 
3 

7 


1 

10 
6 


2 
6 
9 














6 




12 
22 


10 
14 


16 
22 


18 
31 


16 
35 


37 

27 


132 
173 




4 







E 


5 



























6 

7 


I 




















1 




4 


















4 




8 


1 


5 
4 
















6 




9 


1 


5 
2 


1 

6 












11 




10 




1 










9 
























Sum 


7 


17 


26 


24 


41 


25 


38 


49 


51 


64 


342 



Table 4a. 



COS0 





^\ ft 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 


6 
5 


10 
8 


6 

7 


9 

18 


8 
10 


7 
15 


2 

16 


"~13~ 


1 

22 




49 




2 


44 


158 




3 


4 


14 


19 


27 


19 


38 


46 


41 


35 


18 


261 




4 


2 


4 


9 


4 


6 


6 


9 


6 


3 


9 


58 


E 


5 


1 
1 


10 
3 


8 
12 


6 
5 


11 

5 


10 

7 


5 
8 


8 
10 


2 

7 


5 
1 


66 




6 


59 




7 


2 


2 


4 


12 


7 


6 


4 


9 


5 


4 


55 




8 


1 

1 

7 


4 
9 
15 


8 
12 

29 


12 
3 
55 


9 

19 
43 


13 

17 
79 


11 
13 

84 


6 
12 
85 


9 

15 

119 


13 

19 
121 


86 




9 


120 




10 


637 










Sum 


30 


79 


114 


151 


137 


198 


198 


190 


218 


234 


1,549 









Table 4b. 





\v k 
i ^^ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 


6 
4 
3 

1 
1 

1 


9 
6 

8 


5 
2 
13 


7 
14 
21 


6 
9 
9 


5 

13 
25 










38 




2 


12 
27 


12 

28 


21 

21 


43 
3 


136 




3 

4 


158 
1 


E 


5 


5 
2 
2 
2 
5 
2 


3 

7 
3 
2 
1 
2 
















9 




6 


1 
7 
2 
1 
1 














11 




7 


3 
3 

7 
2 












15 




8 


1 


4 
11 
4 










14 




9 


2 

4 








27 




10 




2 






17 
















Sum ._ 


17 


41 


38 


54 


39 


62 


45 


42 


42 


46 


426 









Table 5a. 





^^\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 


6 
12 
10 
3 
4 


24 
20 
22 
4 
8 
4 
5 
5 
4 
18 


26 
23 
37 
16 
9 
9 
8 
17 
18 
27 


21 
19 
41 
16 
13 
7 
9 
10 
12 
47 


18 
11 
41 
9 
20 
20 
11 
14 
12 
56 


12 
15 
38 

6 
10 
13 
18 
14 

7 
86 


11 

17 

34 

5 

10 
10 
16 
10 
27 
87 


7 

26 
31 

8 

4 
12 

7 
21 
22 
90 


1 
35 
10 
2 
4 
4 
8 
8 
25 
137 


1 

23 

10 

7 

3 

2 

3 

15 

12 

116 


127 




2 


201 




3 


274 




4 


76 


E 


5 - 


85 




6 


81 




7 


1 
5 
2 
3 


86 




8 


119 




9 


141 




10 


667 










Sum 


46 


114 


190 


195 


212 


219 


227 


228 


234 


192 


1,857 









136 



Table 5b. 





^\ ft 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 

2 

3 

4 

5 

6 


6 
12 

8 
3 

2 


21 
17 
16 


23 

13 

27 
13 
4 
1 
1 
6 
6 
1 


20 

14 
22 
11 
6 
3 
3 
2 
3 
4 


in 
10 

25 


9 
9 

27 


9 
11 
18 


5 
25 
17 






109 




35 


21 


167 
160 








27 


E 


2 


12 
10 
6 
6 

2 

1 












26 




6 
9 

4 

2 
1 










20 




7 




3 
2 
3 

1 


5 
6 
15 
3 








27 




8 


2 
2 


5 
10 

9 






33 




9 


1 
4 




47 




10 




24 














Sum 


35 


65 


95 


88 


88 


67 


67 


71 


43 


21 


OK) 









Table 6a. 



Sum 



l._._ 
2____ 
3____ 
4-._ 
5__._ 
6._._ 
7„._ 
8._-_ 
9____ 
10... 

Sum 



111 
45 
47 
11 
11 



3 


458 


- 


322 


3 


319 


4 


120 


l 


109 


1 


96 


1 


120 


4 


121 


17 


L46 


85 


490 



237 



256 



228 



183 



2,301 



Table 6b. 





I\! 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




i._. 


107 
41 
44 
9 
9 
2 
3 
7 
6 
4 


80 
30 
46 
11 
8 
4 
7 
3 
5 
3 


58 
39 
47 
12 
13 
8 
9 
6 
7 
2 


39 
52 
37 
8 
8 
4 
5 
9 
6 
4 


48 
17 

48 
18 
8 
5 
6 
7 
8 
5 


39 

24 

29 

21 

10 

6 

10 

7 

4 

4 


29 
34 
1 
2 
20 
12 
8 
5 
8 
8 


18 
28 


14 

21 


3 


435 




2... 


293 




3 


252 




4 








81 


E 


5 

6 








76 




15 
20 
11 
16 
5 






66 




7 


3 
11 
19 
12 




71 




8 

9 


66 




11 

12 


90 




10 


59 










Sum _._ 


232 


197 


201 


172 


170 


154 


127 


113 


80 


33 


1,479 







Table 7a. 





\ k 
i ^\ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1... 



























2 


1 


2 


2 


6 


4 


8 
2 


8 

1 


7 
3 
2 






38 




3 


13 


14 


33 




4 










2 


4 


E 


5 . 








1 




2 

1 
1 




2 


5 




6 




1 






3 






6 




7 








1 

1 

1 

19 


1 
1 
3 
13 




4 




8 








1 

2 
7 






3 




9 










1 

7 






7 




10 _ 




5 


1 


4 


13 


14 


83 












Sum 


1 


9 


3 


12 


16 


21 


26 


34 


31 


30 


183 









137 



Table 7b. 





\ fc 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 


1 


2 


2 


8 
1 


5 


11 
3 

1 
2 
3 
2 
1 
2 
18 


14 
1 


9 
6 
4 

1 






52 




3 


19 
1 

1 
1 
1 
4 
4 
38 


24 

1 
2 
1 


54 




4.__ 








3 


10 


E 


5 








1 

2 


3 

2 

1 

1 

3 

29 


10 




6 _ 




1 

1 




1 
1 
2 
3 
13 


11 




7 . 






2 
3 

5 
42 


8 




8 _ 








2 
3 

48 


13 




9 __ 






1 
6 




21 




10 


1 


8 


10 


213 










Sum .. . 


2 


12 


9 


22 


28 


43 


54 


72 


69 


81 


392 









Table 7c. 





\ fc 
i \. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 


1 


2 


2 


8 
1 
1 

1 
2 


7 


12 
4 
1 
2 
3 
3 
2 
2 

24 


14 
2 


11 

6 
5 
2 
1 
2 
3 
6 
51 






57 




3___ 


20 
1 
1 
3 
1 
5 
5 

51 


27 
2 

3 
2 


60 




4___ 








3 


13 


E 


5___ 








3 

2 
1 

1 

4 

41 


12 




6__. 




1 
1 




1 

1 
2 
3 
15 


15 




7 






9 




8 






1 
1 

12 


3 

6 
64 


17 




9 






1 
9 


28 




10 


1 


9 


277 










Sum 


2 


13 


12 


27 


32 


53 


68 


87 


87 


107 


4S8 









Table 8a. 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 

3 


2 


5 

1 


5 


7 
2 
2 


13 
1 

1 
1 
2 


18 
3 


15 
1 


10 
5 
2 
3 
1 
3 
1 
6 

50 






75 




18 
2 
3 
5 
3 
5 
7 

48 


22 
2 
3 
2 
4 
3 
5 

59 


53 




4 




2 


11 


E 


5 _ 






3 
2 

5 
2 
3 

27 


1 
5 
3 

1 
8 
49 


14 




6 _ 




1 






18 




7 _ 




1 




19 




8 






2 
2 

17 


3 
2 
23 


17 




9 






1 
11 


34 




10 




10 


294 












Sum . __ 


2 


17 


20 


32 


46 


63 


83 


81 


91 


100 


535 









Table 8b. 





\ k 
i v n Ss 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 — 




2 


6 

1 


4 
1 
2 


11 

1 
1 
2 
4 
1 
3 
3 
22 


10 

1 
2 
2 
1 
1 
3 
2 
30 


16 
2 
2 
3 
1 
2 
1 
4 

49 


9 
6 
2 
4 
4 
1 
1 
3 
41 




1 
22 
3 
1 
1 
3 
6 
9 
66 


59 




3 . 




11 
1 

1 
4 
4 
3 
5 
53 


45 




4_ . 


1 




14 


E 


5- . 






13 




6— 






1 


1 


17 




7 






12 




8 






1 

1 

10 


1 
2 
16 


19 




9 




1 
5 


30 




10 


1 


293 










Sum 


2 


8 


20 


27 


48 


52 


80 


71 


82 


112 


502 









138 



Table 8c. 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2... 




2 


5 


4 


11 


10 


16 


8 
3 






56 




3-- 




8 


17 


28 




4... 



















E 


5-_. 



























6 



























7 



























8. . 



























9--. 






















(1 




10 





















































Sum 





2 


5 


4 


11 


10 


16 


11 


8 


17 


SI 



Table 9a. 





i \ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 

2 


1 


1 

2 

1 


2 
5 

1 


1 

4 
2 
2 
2 
3 


2 

8 
6 

5 
1 

1 
3 
7 
9 
3(3 












7 




4 
4 


6 
2 


7 

15 

3 

..... 

3 

4 

10 

50 


8 
10 

1 
4 
4 
6 
9 
6 
82 


11 

( .) 
3 
2 
2 
2 
2 
11 
78 


55 




3 




50 




4 

5 

6 

7 


1 
2 


15 


/•: 


3 


1 
2 
1 

9 
5 
17 


4 
4 
1 
5 

8 
38 


2 
4 
4 
3 

7 
45 


20 
25 




2 

4 
2 

8 


22 




s 


1 
1 

4 


3 

7 
20 


47 




9 

10. 


66 

390 










Sum . 


11 


23 


43 


50 


78 


08 


73 


101 


130 


120 


097 









Table 9b. 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 

2 


1 


1 


2 
3 


1 
2 
1 


2 
6 

2 












7 




2 
3 


5 

1 


6 

10 


8 

7 


11 
6 


43 




3-_. 






30 




4__- 











E 


5 . 



























6 

7 -- 


1 




















1 




1 

1 


















1 




8. 


1 


6 
3 
















8 




9. . 


2 
2 


1 

6 












6 




10 -_ 




1 












g 






















Sum . . . 


3 


4 


14 


8 


17 


5 


6 


16 


15 


17 


105 









Table 10a. 



COS0 



Sum 



Sum. 



46 



73 



3 


1 


4 


8 


11 


12 


1 




3 


4 


5 


5 


4 


5 


10 


9 


9 


10 


70 


77 



33 
62 
58 
23 
37 
46 
43 
80 



120 



131 



139 



Table 10b. 



cos9 





\ k 
i \ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 


3 
2 
2 
2 


5 
1 


5 
2 

1 


3 
3 
3 


4 
3 


6 
3 
6 


2 
3 
9 


1 
6 

7 






29 




2 


10 

2 


15 
1 


48 




3 


31 




4 






2 


E 


5 


4 
2 
2 
3 
2 
1 


3 

8 
2 
3 
















7 




6 . 


1 
1 

1 


3 

7 
3 
2 
2 














14 




7 _ 


2 
5 
5 
3 












14 




g 


4 

7 
3 










19 




9_-_ 


1 

5 








17 




10 _ 






2 






16 
















Sum 12 


20 


24 


26 


22 


29 


20 


16 


12 


16 


197 



Table 11 



COS9 





^\ k 
i \. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 


5 
10 
4 
2 
1 


9 
13 
12 

2 
4 

11 
5 
5 
2 

11 


12 
10 

20 
12 

4 
7 
9 

17 
9 

28 


14 
10 
9 
22 
10 
7 
8 
10 
14 
46 


10 

7 
11 

6 
22 
22 

7 
19 

7 
42 


10 
9 
11 
4 
6 
13 
19 
18 
14 
66 


5 
5 
6 
2 

10 
10 
17 
31 

74 


6 
12 
9 
2 
3 
5 
9 

20 
29 
92 


5 

18 

6 

1 

3 

4 

5 

11 

14 

119 


"""!§" 

4 
3 
2 
4 
5 

10 
19 
115 


76 




2 


107 




3 


92 




4 


56 


E 


5 


62 




6 


83 




7 




77 




8 


1 
3 
1 


128 




9 


142 




10 


594 










Sum ._ 


27 


74 


128 


150 


153 


170 


167 


187 


186 


175 


1,417 









Table lib. 



COS0 





>V \ k 

i \. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 . 


5 
10 
3 

2 

1 


9 
10 
12 
1 
2 
1 
3 
1 
2 
3 


11 
3 

18 
10 
4 
2 
1 
6 
4 
1 


13 
3 
5 

18 
5 
5 
3 
3 
4 
2 


9 
6 

8 


7 
3 

7 


5 
3 
2 


6 

12 
3 


4 
16 




69 




2 


12 


78 




3 .. 


58 




4 . 






31 


E 


5 


13 
16 
5 

7 
1 
2 












25 




6 


8 
11 
4 
4 
3 










32 




7 




6 

10 
13 

4 








29 




8 


1 
1 


5 
7 
12 






37 




9 


1 

7 




37 




10 




34 












Sum 


23 


44 


60 


61 


67 


47 


43 


45 


28 


12 


430 









Table 12a. 





\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 


81 
38 
19 
17 
16 
7 
7 
11 
13 
4 


62 
30 
23 
20 
20 
12 
17 
9 
10 
19 


50 
38 
22 
19 
19 
15 
19 
13 
16 
23 


26 
57 
30 
16 
22 
12 
12 
19 
11 
38 


29 
19 
46 
26 
17 
10 
15 
21 
20 
44 


25 
13 
18 
38 
20 
14 
22 
18 
11 
52 


10 
16 
3 
7 
29 
33 
17 
19 
23 
73 


15 
16 
3 

25 
36 
28 
28 
71 


10 

14 

4 

2 

1 

5 

7 

20 

38 

93 


8 
11 
3 


316 




2 


252 




3 

4 

5 

6 


171 
145 


E 


1 
2 

2 

18 
99 


147 
135 




7 

8 


154 
160 




9 

10 


188 
516 




Sum 


213 


222 


234 


243 


247 


231 


230 


224 


194 


146 


2,184 
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Table 12b. 







1 


2 


3 


4 


5 


G 


7 


8 


9 


10 


Sum 




1 


77 
34 
18 
15 
14 
6 
5 
9 
11 
2 


60 

27 
20 
17 
17 
9 
10 
3 
8 
6 


50 
35 
20 
15 
16 
10 
15 

10 

4 


26 

53 
29 
15 
14 
6 
9 
8 
(i 
8 


27 
18 
45 
25 
11 
7 
10 
10 
9 
6 


25 
8 
17 
37 
15 
9 
14 
10 
6 
7 


10 

13 

1 

5 
29 
L\ r > 
8 
9 

14 
12 


15 
14 


10 

J 3 


8 
10 


308 




2 


225 




3... 


150 




4... 








129 


E 


5-.. 








116 




6 

7 


21 
31 
19 
20 
9 






93 




5 
17 
28 
20 




107 




8 




93 




9 


11 
21 


123 




10.-. 


95 










Sum 


191 


177 


183 


174 


168 


148 


126 


129 


93 


50 


1,439 









Table 13. 





^\ k 
i ^^ 


1 


2 


3 


4 


5 


6 


i 


8 


9 


10 


Sum 




1 


5 
6 
5 
1 
2 
2 
1 
2 
3 
6 


15 
5 

10 
3 
7 
6 
3 
7 
4 

10 


6 

7 
9 
8 
3 
7 
8 
3 
5 
24 


8 
13 
11 
3 
5 
6 
7 
10 
6 
21 


5 
5 
9 

9 
5 
11 
9 
13 
38 


6 

10 
5 
6 
4 
14 
10 
15 
19 
59 


5 
10 
9 
2 
4 
12 
10 
20 
7 
G2 


2 
6 
6 
3 
2 
8 
9 
10 
18 
57 


1 

12 
5 


3 

21 
10 
1 
4 
4 
7 
4 

13 
106 


56 




2 


95 




3 


79 




4 _. 


33 


E 


5 

6 . 


2 
6 
6 

7 
10 

98 


42 

70 




7 

8— 


72 

87 




9 


98 




10 .- 


481 










Sum 


33 


70 


80 


90 


110 


148 


141 


121 


147 


173 


1,113 









Table 14a. 





i ^^. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 

3 

4 


3 

1 


3 


11 
1 
1 


5 

1 

1 
1 
2 
2 
2 
10 


15 
3 


16 
2 
2 
1 
2 
1 
3 
4 

22 


15 

2 
4 
2 
1 
3 
1 
6 
30 


11 

8 
3 

2 
7 
1 


3 

13 
2 
1 
1 
3 
5 
5 

50 




82 




23 
3 
2 

4 
5 
3 
7 
50 


54 






16 


E 


5 ... -_. 






2 
5 

3 
2 
2 

14 


11 




6 




1 


1 


23 




7 




18 




8 






2 


18 




9 




1 

4 


3 
31 


30 




10 


2 


6 


219 










Sum 


6 


9 


22 


25 


46 


53 


64 


66 


83 


L97 


471 









Table 14b. 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 . 



























2 _ 


4 

1 


3 


12 
1 
1 


2 
2 
3 
10 


16 
3 


18 
4 
2 
1 
2 
1 
3 
4 

22 


16 
2 
4 
2 
3 
3 
1 
6 

31 


17 
10 
3 
2 
7 
3 
1 
3 
33 


3 

16 
2 
1 

2 
3 
5 
5 

51 




95 




3. 


28 
3 
3 

5 
3 

8 
51 


66 




4 




16 


E 


5— 






2 
5 
3 
3 
3 
14 


12 




li 




1 

1 


1 


26 




7 _ 




21 




8.-- 




2 


20 




9 




1 

4 


33 




10 .. 


2 


7 


225 










Sum 


7 


10 


24 


27 


49 


57 


68 


79 


88 


105 


514 









512908—59- 
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Table 15a. 





\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 




2 


2 
2 
1 


3 

3 


4 


2 


5 


3 

1 
2 
3 
2 
1 
3 
3 
23 


2 
6 
2 
1 
4 
3 
3 
9 
28 




23 




3 




21 

4 
2 
1 
1 


33 




4 




1 


1 
2 
2 
1 
1 
2 
17 




2 


13 


E 


5 






8 




6 




1 


3 

1 




3 
3 

7 
3 
26 


_____ 

3 

5 
17 


16 




7 






13 




8 




1 


4 
1 
14 


22 




9 






10 
21 


33 




10 


2 


5 


6 


159 




Sum _. _ __. 


2 


10 


15 


25 


30 


44 


35 


41 


58 


60 


320 









Table 15b. 





^\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 


1 


2 

1 
1 


3 

2 

1 


5 
4 


8 
1 
1 
4 
4 
2 
1 
3 
18 


2 

1 
3 
1 
3 

4 

7 
3 
29 


5 


6 
2 
4 
3 
3 
1 
4 
3 
24 


3 
11 

2 
2 
4 
3 
3 
10 
30 


27 

12 

27 


36 




3 


49 




4 




4 
3 


20 


E 


5 




2 
1 


19 




6 _ 




1 


3 
2 

1 


20 








4 
3 

6 
20 


18 




8 




1 


4 
3 

19 


26 




9 




40 




10 


2 


6 


7 


182 




Sum 


3 


12 


19 


38 


42 


53 


45 


50 


68 


80 


410 



Table 15c. 





i ^\. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 

3_ _ 


3 


4 

1 
1 
1 

1 


4 
2 
1 


6 

5 


8 
1 
1 
4 
4 
2 
1 
3 
15 


2 
2 
3 
1 
3 
4 
3 
2 
30 


5 
1 
4 
3 


4 
5 
3 
3 
2 
4 
2 
21 


4 

14 
2 
3 
4 
3 
3 
5 
28 


1 
33 
3 
3 


44 
63 




4 




20 


E 


5 




3 
2 


21 




6 




3 

2 

1 


20 




7 




4 
1 
4 

18 


2 
2 

7 
24 


19 




8 








15 




9 _ 






2 

14 


25 




10 


2 


7 


6 


165 




Sum 


5 


15 


19 


32 


39 


50 


40 


51 


66 


75 


392 



Table 16a. 





^\ k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 




3 


4 


4 
3 


6 


3 

1 


4 


2 
1 
1 
1 

16 
8 
2 
5 

13 


1 

7 

"T 

14 

8 
4 
4 
25 




27 




3 




21 
2 
1 

7 
11 
5 
6 
22 


33 




4 








3 

1 
9 
5 
5 
5 
15 


6 


E 


5 










1 

10 
4 
2 
2 
14 


1 
10 
6 
1 
5 
15 


6 




6 

7 


2 


1 


3 

1 
1 
1 
2 


8 
1 
2 
3 
11 


80 
44 




8 




1 


23 




9 

10 


1 
1 


32 




3 


121 




Sum 


4 


8 


12 


32 


39 


42 


47 


49 


64 


75 


372 
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Table 16b. 





S \ k 
i ^^^ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2 


1 


6 


10 


7 


16 
3 


10 
1 
2 


13 

2 

_____ 

8 
6 
6 
6 

20 


13 

7 
1 
1 
17 
10 

i 
23 


2 
13 
2 
2 
13 

4 
4 

27 




78 




3.. _ 


13 

1 

""16 

14 
6 
9 

28 


39 




4 

5 


1 






1 


8 


E 


1 

1 






8 




6 _ 


2 


3 
2 
1 


8 
4 
6 
5 

12 


9 
2 
3 
6 
11 


7 
9 
2 
6 
23 


78 






52 




8 






35 




9 


1 
2 




38 




10 


4 


3 


153 










Sum 


7 


12 


19 


43 


50 


60 


65 


s 1 


72 


SI 


is.) 



Table 16c. 





\ fc 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1- . 



























2. _ 


5 


8 

1 
1 
1 


11 
3 


12 
9 


14 

4 
3 
2 
2 
3 
3 
3 
10 


6 
3 
1 
3 
4 
4 
4 
5 
15 


13 

5 
8 
4 
3 
4 
3 
2 
22 


13 
14 
5 
3 

4 
2 
1 
4 
15 


4 

34 
2 
4 

10 
4 
5 
8 

21 


2 
56 
4 
4 

1 
4 
4 
9 
16 


88 




3— 


129 




4__. 


2 


26 


/•: 


5__. 




3 
2 
2 
1 
2 
13 


24 




6 






28 




7... 






23 




8... 






2 
2 
2 


23 




g 


1 

2 


1 


36 




10 


117 




Sum 


10 


14 


20 


44 


44 


45 


64 


61 


92 


100 


494 









Table 16d. 





i ^^\^ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Sum 




1 



























2.. 




3 


4 


4 
3 


6 


3 

1 


4 


2 

1 
1 
1 
2 
1 


1 

7 




27 




3. . 




21 
2 
1 
1 
2 
1 
3 

11 


33 




4 










3 

1 

1 


6 


E 


5 










1 

1 
1 
2 

1 
9 


3 
2 

1 
3 
11 


1 

3 

2 


6 




6 




1 




1 


13 




7 






8 




8 . 








1 


1 

1 
14 


6 




9 


1 

1 






1 
8 


4 
15 


14 




10 


1 


1 


8 


79 










Sum 


2 


5 


5 


17 


21 


25 


25 


17 


33 


42 


192 
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